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ABSTRACT 


The disastrous results of large rock slides have 
resulted in the geotechnical world's increased desire to 
predict the degree of stability of rock slopes. The problem 
the engineer faces in making such predictions involves the 
reliabidityyotythe stabidity analysicid=Inethestanalysisheot 
rock slopes, safety factors based on laboratory shear tests 
may, in some cases, be in error by as much as 700%. 
However, in a geological setting like the Frank Slide this 


error may only be in the order of 15%. 


A relatively good correlation is obtainable between a 
mathematical characterization of the microscopic roughness 
of artifically prepared surfaces and the friction angle, Q. 
The extension of this concept to the field scale is hindered 
by the limited number of suitable surface exposures 


‘available for the measurement of roughness. 


The frictional resistance after large shearing 
displacements on natural rock surfaces depends on the 
initial surface roughness and therefore is not a unique 
parameter aoe thie residual strength iS hole Giays: 


Furthermore, the shear strength on _ rock surfaces is 
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independent of sample size if the degree of roughness does 
not change with the area. The field shear strength of 


flexural-slip surfaces may not be at their ultimate value. 


The geological setting of Turtle Mountain differs 
Significantly from what has been shown in many publications. 
The tstLuctunad leform iofrthis *mountain iis can tan ticbhine fenot Tha 


monocline. 


Back-analyses of the Frank Slide indicate that the $ 
requived forharsatety ssactoniofd1i0 isirtsimular tothe Thasic 
fprction angie; O.- The use of peak strength parameters 
from flexural-slip surfaces in the back-analysis results in 


a@tsatety “factoriof 0.86. 
The mining at the base of Turtle Mountain decreased the 


safety factor of the sliding mass by 1% as shown by a Finite 


Element stress analysis. 


(v) 


ee 
sano teod Tang yasn at nvods) desd mbit! sadw work yLinaDeRaple 
& tom venalbtens ite = stusaeda eth 30 10 Renee ate 
; nd 7” 
A odd dad) otsDihab obife anast ods To epeyieie-soea’ 
atest odd of vefinke 2¥.0,1 20 rodbak ydatae 5 aod bextuper 


Syetemezeq dtproite Assq 20 sou ed? = -,h ,otpaR not soba? 


ot etivest aievlens-xosd ent of ads true qile- Iswwxet? moat 
bs: ‘488-0 20° ro#o8% tae 


gid Peanetosh eistavot afew? Iolersd edt 46 pnieim edt 
sditetd wo yd awode as Pt ye aesm vAERite aH% 20 ro2902 ysetee 
nieyisds eeerte! thebeats 
a ete 


ey ; 7a) 7 1 mina > ‘1 : 7 e a 


= 
7 
" 


q 


ACKNOWLEDGEMENTS 


THEE author #sorlie@thas Wthesisemiseqrateful to Professor 
N. R. Morgenstern under whose supervision the research for 
this thesis was conducted. His topic suggestion and 


continued encouragement have been deeply appreciated. 


The author is also indebted to Professor D. M. Cruden 
for his leadership on the Turtle Mountain geological field 
trips and for his helpful guidance in the interpretation of 


the related data. 


The assistance with the work on Turtle Mountain 


provided by D. K. Norris is appreciated. 


Appreciation is extended to Rk. Gitzel under whose 


direction the optical roughness plotter was built. 


The editorial assistance rendered by W. D. Roggensack 


is acknowledged with appreciation. 


The financial assistance provided by the National 
Research Council of Canada, the Department of Energy, Mines 


and Resources, and the University of Alberta is gratefully 


(vi) 


sna 'aanonatpehs Seok as ecla  e 
sHorsionrge vgeeb aged a¥nd Sone 

iby? «HM 1d TokestoTI OF en eels ab Pere 
Biek? Laoipofooep aistavoM ofti0T ode no @tdatebset etd . 
te hoktavesysetnt ed? at enaebiite tabetha eka sot ban watt 
ined oosphan ea : 


Atedavew o4tdut noi Atom ody div somayaisee: ot ve 
nd | | betetosacgs ef @H 0m Matt ‘ae 

. oad 

ea0nw tehao. Lastio 2 oF hebmetxe ef aeitsisotqga J a 
thiud eau Chait eeandpiot Leottqo oft etal) 


“ 


‘einen a e yd boveboot constetees Lsitorthe edt 


acknowledged. 


Finally, the author wishes to acknowledge, with many 
thanks, the continual encouragement, support and editorial 
assistance rendered by his wife, Barbara Ann, as without 


her, this thesis would not have been possible. 


(vii) 


*% Y anu 
: ‘ t - . : 
Se iy 
7 =~ 7 
os ie ay ' 
ie eed 5 “fa 
ips er an A 


TABLE OF CONTENTS 


Release Form 
Title Page 
Approval Sheet 
Abstract 
Acknowledgements 
Table of Contents 
List of Tables 


List of. Figures 


GHA TE Ra INTRODUCTION 
as, General 


WZ Purpose and scope 


GHA Pwr Ree THE STABILITY OF ROCK SLOPES 
ZN introduction 


Dae Analytical slope stability 
analysis 


2a3 Physical scale models 


2514 Governing parameters in slope 
stability analyses 


2a Case histories 


2.6 The reliability of rock slope 
stability analyses 


(viii) 


be 


Ve 


i a 
oh | - rent ose 
wi | pa 
iy . ationspbelwontoe, 
(ritv) 7 | a ainetned 30 older 
(bee) a? sefdsT to veht 

(thixy eemmpit To fete ok 


HOT TIUCOATUT ‘a al 


i”, A | | ' fenensd T.T 


eo egoot bas stogust %.F 7 


t aia) 
exioue AD0A AO WTaTeaTe TdT Ti AeTSAND iP 

+ nae 

ee | soitorbordnr f° .S8 , = 


qabtiaate aqola Leot> ma) See ee 7 
; at £1 6 bh aa 


as (salaet isole fe Le) | 


CHBPIRR I 


CHAPTER LV 


TABLE OF CONTENTS (CONTINUED) 


THE SHEARING RESISTANCE OF NATURAL AND 
ARPTIFICHAL ROCK SURFACES 


Br 4 


35 2 


Introduction 

The influence of surface 
irregularities on the shear 
strength 


The effective i on natural rock 
surfaces 


Basic? trtetion@ angle’ of rock 

The effect of microscopic 
roughness on the basic friction 
angle 

Definition of roughness 

Sample preparation 

Shear tests and results 
Digitization of roughness profiles 


Digital computations 


Amplitude distributions of 
asperities 


Roughness characterization versus 
frictional resistance 


Duscussion of results 


Page 


oid. 


SHEAR TESTS ON NATURAL ROCK DISCONTINUITIES 
FROM TURTLE MOUNTAIN 


Introduction 
Sampling location 
Sampling procedure 


Sample preparation 


(ix) 


82 


OD 


83 


86 


Phe hols Lexvted no 4 Anau = Bi iy Ae erg 


A501 to efpms wolto fat ofesh aye 


rae | > on 
no os, ge _ |<? i 
eeondguey to noftiniteG Of : ia 

noktsrsqszq okqmse v.£ F “© aS v4 

ativass has aves ‘teeta (Be ee en 

aelitorg a2endpuvor to nokissitipig  C.€ ‘ - | 
xpoibhedkeinn's fetinsa Oat 7 

+0 See) ree - iS 

aveisv dank a a 7 St .é : 7 | 
sfee0% Ieaoks Ae 

“as toe2 Ie aananeaeea Er.e : ~ ‘yt a | 

Pay 


-eaETTOMT mode i AOO% TK 7 


CHAPTER Y 


CHAPTER VE 


CHAPTER Vit 


TABLE OF CONTENTS (CONTINUED) 


4.5 Testing apparatus 


4.6 Testing program 
Ah pe, Results of shear tests 
8 Discussion of results 


THE GEOLOGY OF TURTLE MOUNTAIN 


eal Entro duction 
B02 History of investigation 
Be Succession of rock units within 


the Turtle Mountain area 
5 4 Geology of Turtle Mountain 


Sie: Discussion 


THE MECHANICS OF THE FRANK SLIDE 
6.) Inecoduct ion 
oi 4 The stability analysis 


653 The effect of mining on the 
stability of Turtle Mountain 


6.4 Discussion 


CONCLUDING REMARKS 


LIST OF REFERENCES 


APPENDIX A 


OPTICAL ROUGHNESS PLOTTER 
Neca Introduction 
Aw Theory of operation 


A. 3 Construction 


(x) 


Page 


120 


122 


ast 


oe peice 
pec ibe wa's 


atesavem er to ypotosd 


‘ , ee 

Jarae WKAR SRT IO 2IDWARDMM SET 
7. nobtowboatel = t.8 
--afeyisnn ytitidete oft S48 


APPENDIX B 


APPENDIX C 


TABLE OF CONTENTS (CONTINUED) 


A.4 Operating procedure 
ne Results of trial tests 


A. 6 Performance 


SHEAR TEST DATA 


SPECTRAL DENSITY COMPUTER PROGRAM 


(xt) 


Page 
186 
arom 


188 


196 


230 


yi MY 


: - : ; te My 7 
ye | tt BAY vs ‘i i mn : if ; 
es a : \. 


nabs we i 7 + 


be 


Lietii: =: Af Gin 
4 a 


List OF PABLES 


Summary of case history studies 


Sample designation for roughness study 


Roughness parameters for surface profiles 


Summary of shear strength parameters 


Width of cracks on Turtle Mountain 


Results of stability analyses on Frank Slide 


Shear test data 


(xii) 


Page 


56 


100 


130 


160 


197 


oat 
| ; ter 


 ebti? {184% 0 enoyions) ee 20. attveed ~ 


 ptpb tees Tye 


Figure 


Lio. -OF FIGURES 


Limit equilibrium analysis of a sliding block 
Livwstratzon of circular are fallure surface 
Illustration of a general failure surface 

A three-dimensional wedge failure 
PlLustration of @ potential failure surface 
which includes pre-existing joint sets and 


Tatacterocek 


Failure envelope for specimens with flat 
surfaces 


Failure envelope for multiple inclined 
surfaces 


An example of a natural discontinuity 
illustrating first- and second-order 
irregularities (after Patton, 1966) 

An example of a time-history record x(t) 


Casting of limestone samples 


A typical step loaded, shear load versus 
deformation plot 


Failure envelopes for rock surfaces with 
various roughnesses 


An example of a surface trace from the 
Talysurf 


Surface profiles for the dry sanded sample 


Surface profiles for the diamond-saw cut 
surfaces 


Surface profiles for sample lapped with #45780 


iat 


(xiii) 


Page 


24 


29 


25 


26 


26 


emt) 


58 


a9 


59 


60 


61 


62 


63 


64 


65 


66 


as 


ve 


“a 


* 


woold piibite « Yo siaiitions auttarctape sek 
ede te aaa ous astepato 2% nottertew Lily $.¢ 
onstine aru fked Estanep 5 te polgexdao Ele £.$ 
etl te? sales se See he 


sxeitne exvtist ae . a PET 
bas 2tea totot pais . Res’ 


tart dtiw eremtoege 120% egolovie at ul iad a | 


a 


fearlouvt ofq is Lue 3103 seolevne eqerss ~~ SE 
e2aneiwe . 


ytionbemovech Lsaotba. 5, to Sfate<e mk €.6 ate 
, awebte-bhowes bas -FenFd, py ee a 
(@0er ,nodtet aotts) 2 > 
(#) x fxosst ytoteld-smit 6 to etme wh 6) oie 
| q 

asiquse saotasmil to anisteso o.6% 
svezev fsol raede ,bebsol qoete Lsot Ce : 
70 J not as if 
dtiw 2e0stive AQO7 10> Seendovo3 etut tet ¥.€ 
asnusnipuot eeortsv > w 
adt mot? sasts spstme & 20 siqmsxe ms §.€ 


eiqnes bebmse y25 sat 20% aol: tory oosta0e 
too wneebabebth Sit/.29% mel ttoag: ooebaee 


OR\28# Ativ boqgst Sigaea tot aaterorg soeaave 


LESTOOF FEIGURES ((CONTENUED) 


Page 
Surface profiles for sample lapped with #200 
Clit 67 
Surface profiles for sample lapped with #400 
Grit 68 
Surface profiles for sample lapped with #600 
Girt 69 
Amplitude distributions of asperities for dry 
sanded sample 70 
Amplitude distributions of asperities for 
diamond-saw cut surface qe 
Amplitude distributions of asperities for 
surface lapped with #45780 grit V2 
Amplitude distribution of asperities for 
surfaces lapped with #200 grit ve 
Amplitude distribution of asperities for 
surfaces lapped with #400 grit 74 
Amplitude distribution of asperities for 
surfaces lapped with #600 grit a 
Centre-line-average (CLA) versus the friction 
angle 76 
Roughness parameter, Zo, versus the friction 
angle 148) 
Roughness parameter, Za, versus thestr3 ction 
angle 78 
A typical autocorrelation plot FAS 
A typical spectral density plot 80 
Limiting spectral criteria for roadways (after 
Houmbowe, 15617) 81 
The sample block within the slide debris 101 
A flexural-slip surface within the sample 
block 101 
Coring on the sample block 102 


(xiv) 


- _e | > 


4 


ef 


as 


— - 
> - 
7 7 - rf = 


eis 2d sotstroans to | ar t 
202: asisttonas to ‘estan | ar.€. 
wa 7 { 7 
2 Bene tr.t 
a? ake 29. aos gue aga south 


+E 54 art 
ser Ee a aad spe, 


t0% maine = ae a ae ampere ~ ere . iq 


at ee ta ssid oT |, | os .é 


RoLpotazt sd+ ataqv. Cea d5y) ame <8 


fofiots% oft auetev. oS vtefemsrag a ie ESE 


- toLq act+eferzosetue Inpkay?. a 
tog) vets, isti>oqe ptank ‘ 


nottaList sdt enatov y,8 rere ery wage ¢g.€) 


Liat OF FIGURES®= (GONTINUED) 


Figure Page 
u.4 A ten-inch diameter cored sample 102 
4.5 The modified shear box for the two-inch 

samples 102 
4.6 The Wykeham-Farrance shear box 103 
uh Natural discontinuities 104 
4.8 The small roughness measuring device 105 
4.9 A joint surface after shearing 105 
4.10 A typical shear load versus deformation plot 106 
4.11 A typical vertical deformation versus 

horizontal displacement plot 107 
4.12 Failure envelopes for the bedding planes 108 
4,13 Failure envelopes for the flexural-slip 

surfaces 103 
4.14 Failure envelopes for the joint surfaces 110 
4.15 Failure envelope for surfaces lapped with 

#45780 grit ns 
4.16 Failure envelope for the diamond-saw cut 

surfaces Taiz 
4.17 Natural surface profiles of a bedding plane 

(sample B-3) 113 
4.18 Natural surface profiles of a flexural-slip 

surface (sample S-2) 114 
4.19 An idealized surface showing the effects of 

different sizes of specimens (after Deere, et 

ae, oo) vale 
el Vocation of Turtle Mountain 13a 
eee Cross-section through the north peak of Turtle 

Mountain (after Daly, et als, 1912) 132 
De 3 Allan's cross-section through the south peak 

Of Turtle Mountain (ater Aldan, 1933) 133 


(xv) 


oo.tve® sieinesactals penance 
paiawede goats yatis setot & 


tol nodten seb everaw head shite fastqy? A 


2ue20v aotten oteb L: 
toly tema: 


#onnla piibbed ox? 702 Beqofovee ordi ies 
giia-Larvxelt out to? eegelevas at 


Mt tf 


esvaetiue snot air 108 ‘deudlinwis jodiibes 


sh pte athe 10% vee a 
tum ere aoy Bat Ot dqotexce 3: 
aAD: 


wireta putbisd s 0 gefitord eoetaue a 


4 


| ahi tantiness 8 to aebitorq eosk hcaaae 
ee 2 of =e) 2 
aries woke sansa? 


ar) 
ay 
oro 
the 


sre 
: er.a 


ar a 


aren - 
ate 
tr ~ 


es 


Figure 


LIST OF FIGURES (CONTINUED) 


MacKay's cross-section through the south peak 
of Turtle Mountain (after MacKay, 1932) 


Norris' cross-section through the south peak 
of Turtle Mountain (after Norris, 1955) 


Airphoto of the Frank Slide area 


A plan showing the general areas of the 
geological traverses run 


Views,ot north face of Hillcrest Mountain from 
Turtle Mountain showing the Turtle Mountain 
anticline in profile 


The crest of Turtle Mountain 


The Banff Formation in the "hoodoo" weathered 
area 


Oblique aerial view of the Frank Slide fron 
the northeast 


Striations on bedding surfaces in an abandoned 
limestone quarry at the north end of Turtle 
Mountain 


A rock within the slide debris showing a 
sheared surface 


Geological map of Turtle Mountain area 


Cross-section through Turtle Mountain along 
line A-A' 


Cross-section through Turtle Mountain along 
Laine, B= BS 


Cross-section through Turtle Mountain along 
Lane C-C! 


The average normal and shear stresses in 
relation to the failure envelopes 


The effect of cohesion on the factor of safety 
LOG Section B=3* 


The Finite Element grid for the stress 
analysis of Turtle Mountain 


(xvi) 


Page 


a Rela) 


138 


138 


139 


140 


144 


145 


she" 


sys: 


mort ghetmoMm tapipls 
he wtinon oftitt mb 


eet 
a nied nv? al taut tp sient (Bab 
 . feredktsew Pouboad™ ant of nbséaz07 vast oi ‘ae or.2 
cE aa. 
fort obil2e amend ott) to sotv Lehr iat hee! 
oat . weak ree eT} a 
benobnsds 15 ck sepstave patbhed ao pact how gree 
gitto? to bre dttom sit Hs ees td S09 Th. a 
rar - (PLR TAVOR fu '/- aan 
| p patwods etideh efite sft ainziW oer A ff 2, 
ter oyetite feteete on 
Sat sous abs*iueM gitar? to gen isodpotoa® whee ie 
eh | profs iieoiua ofaa dpoo mds sot hpea-2a0%9 ar. 


phols Stam ott? dpb wt siete 
Her | - Amo suit 


Figure 


6.4 


LIST OF FIGURES (CONTINUED) 


The initial stress distributions along the 
slip surface 


The Finite Element grid around the mine area 
with the mine open 


The changes in stress distributions 
The inward movement of the mine walls 


The style of deformation of Turtle Mountain 
due to the mining 


Schematic diagram of light probe for the field 
roughness plotter 


Schematic diagram of holding frame for the 
field roughness plotter 


The optical roughness plotter 


Control panel for the optical roughness 
plotter 


Roughness plotter set-up in the field 


Profile tof artificialisurtadce shown, in lLefgure 
his 


Roughness profiles from the field set-up shown 
LUPE Ugure yA. 5S 
Shear load and vertical deformation versus the 


horizontal deformation 


Cosine taper data window (after Bendat and 
Piersol, 1971) 


(xvii) 


Page 


194 


nos 


195 


23.5 


eer 


EFS 


‘aga tans = pers edt 
ailsv one ait to sneneven, brewer oat 


prea we 9IF10T te mo aden 20% 


bleft eat sz0t guorg ees ia cae, 


aetsoig veo 


eth vitsmedre 


ght xokY amex? paibiod ta) as 
Ssndpgod Slots 


19220Tq 2 
taste lq aeentpior sniaehds iia 


azonodpyos Lspiigo sgt tot fsnsq for sod 
aettola 


bLert od? at qu-$82 rotioly axen dpro8 


ome £7 it fwoda epsttue Istoriistrs to otizoxs 
«4 


avons gu-tee fleii ons moth nel maas eaZond et. 


-4 siued 


nodose aot 


| odt aulatov notsemt0teb fsoistiev auf ieee bias 


bas +sbeod 1925) vob esr 1eded 2 1ba0D 


Ay ver 


“Bo Ee. 


a 


CHAP |S 


INTRODUCTION 


1.1 GENERAL 

MVEELAGESSCRUMBRE INOSETALIAN® ELOOD." This was the 
headline of an article which appeared in the Edmonton 
soutnal onfPOctobert10;2196a;eand #ircrchatedst he) disasters -to 
which the villages along the Piave River, Italy, had 
succumbed. A gigantic rock slide just upstream from the 
Vaiont Dam plunged into the reservoir generating a wave of 
water 300 feet above the crest of the dam. The wave toppled 
over the dam crest flooding the Piave River valley. As the 
flood swept down the valley, several villages were destroyed 
and approximately 2,609 lives were lost (Kiersch, 1964; 
Muller, 1964, 1968: Mencl, 1966; Skempton, 1966; Kenney, 


19672" Nonveiller> 3967) 


Throughout history, similar events have occurred all 
over the world. To mention a few, there was the Flim Slide 
in Switzerland Ol(Heim .e4932)4f the Tatjord Slade in Nonway 
(Bugge, 1937), the Madison Canyon Slide in Montana (Witkind, 
etal, A962)" and the: Frank SdidesamaAdbertasd(Daly, eb eal., 
Poa) In all these cases, numerous lives were lost which 
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potential slide areas. 


As nations throughout the world are becoming more 
industrialized, there is an ever increasing need for more 
natural resources. Japan, for example, has in recent years 
bought coal from Alberta and British Columbia making it 
economically feasible to mine this product in large 
quantities. In some areas of the world, certain minerals 
have become a scarce commodity and as a result, lower 
guality ore bcdies are becoming more economical prospects. 
In the case of open pit mines, the necessary removal of 
large quantities of material has become economically 
possible today with the advent of large-volume earthmoving 
equipment. The consequence of all these developments is 
that open pit mines are being planned and operated to depths 
and sizes that had not previously been considered (Stewart 


and Kennedy, 1971). 


In these large open pits, one of the more important 
factors which must be considered in an economic analysis is 
the maximum slope angle at which the pit can be excavated 
(Steffen, et al., 1970). This is quite understandable since 
the mintmizatkion of jthe amount of waste rock, which ‘has to 
be excavated in recovering an ore body, requires that the 
ultimate slopes be cut to the steepest possible angle. 


Oversteepening, on the other hand, results in slope failures 
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which are not only dangerous to working personnel but may 


also cover up valuable ore in the event of a slide. 


For many large pits to be profitable, the slopes must 
be designed at an optimum angle. Brawner (1971a) points out 
that on three recent projects with which he was involved, 
each degree of slope was associated with a potential cost or 
saving of from five to 15 million dollars. TALS 
demonstrates how very critical the role of the slope angle 


can be in the economic operation of the project. 


Stewart and Kennedy (1971) have also shown that it is 
not only the ultimate slope angle that has an effect on the 
overall profitability. On the basis of cash flow thev_ show 
that it is often economically more advantageous to use steep 
slopes *iduring ithe Sinitial stripping stagiesah This again 
points to the fact that the mine profitability is influenced 


by slope design considerations. 


1.2 PURPOSE AND SCOPE 

The above discussion indicates the need and 
desirability of accurately determining the degree of rock 
slope stability. The problem the Geotechnical Engineer is 
facedaawiths 2neefufiblingwthus needins whethenton noe ates 
possible, on the basis of present day analytical methods, to 


predictethisadegreesofystabilityato athetydesiredwelevel sof 
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confidence. 


Tte iso there purposelteof oithislethesisir to examine! this 
problem by firstly inspecting numerous documented records of 
failed slopes (Chapter 2) and secondly, by performing a 
detailed back-analysis of the Frank Slide which is located 


in southwestern Alberta. 


The reliability of stablity anaylses depends to a large 
extent on the facility with which the operative field shear 
strength can be established. Various methods have been 
proposed for determining the shear strength; however, no one 
particular method has been generally accepted for routine 
use. These methods are reviewed in Chapter 2 and are 
further investigated by a study of the relationship between 
the microscopic surface roughness and the shearing 


resistance. 


Obtaining a proper case history study requires an 
accurate knowledge of three major items, namely: (1) the 
geometrical and geological setting, (2) the shear strength 
parameters along the slip surface, and (3) the water 


pressures acting within the slope. 


In more detail and with special reference to _ rock 
Slopes, the geological setting involves establishing the 


lithology and stratigraphy, the attitude, geometry and 


mat 


We dens akeod? ats ay. talk “ik yh #0. 
to ebr0de7 ‘Seodnomdo) epodoqun pakersqait yLtertt vi wotsore = 
® patutotrey Yd «yibroves Bas (8 Tetaed9) seqose berks? 
betavel at doLdv obife dasxvt odd to. efeylens-doad Bbeliezed. 
| eeeed hh atetaowdsooe nt 


4 J 
4 


spin 6 oft abiaqeb soalyans volidete Yo etitidsiies ) ) 
vs6de biett evite7eqo 64% dozdw ide ie yriitos2 eds so tne 7x8 
aded sve -a2folgjen enotusy .botekidstse ed nao Aspoe78 
no 6A ,xevavod ;itpnetie tsene od? priniatstob wot herocotg 
saivvor sot baiqesos tilsteasp osed esd bostoe aelevitisq 
eig baa = zetqgsis nt beweiver e1s ebotten seed? od 
neewtad. tiene iisiet edt to ybude # yo betspltesval teddaet 
pritsede? o4) fons sesntippor  sontave staovactoie edd ~ 

some tatee? 

is fy 

as #s67vtvpes yBpte yaotetd a3a259 Iseqotq & pataisido 

edt (PM) tyfonrn acett 4otse ssid? to spbeiwoat etezwOSs 

Adonorte teede odt (8) \paRit+ee Lentpofeep Bre Lentadeaoep 

tetew oft (€) Baw senettue gifs eit pucks arsteasisd 
oqoke odd nidtiw pattos cervesemg | 


Spatial distribution of discontinuities and how these relate 
40 the fatlure surface. The geometric setting involves 
ascertaining the surface profile prior to failure. 
Fstablishing the shear strength parameters requires that 
either field or laboratory shear tests be performed. A 
knowledge of the water pressures within the discontinuities 
is important because of its relation to the effective 


stress. 


With these requirements in mind, the investigation of 
the Frank Slide is organized and presented as follows: 

(1) Chapter 4 is devoted primarily to the determination 
of the shear strength parameters used in the back- 
analysis of the Frank Slide. 

(2) Chapter 5 gives the geology of Turtle Mountain and 
the geological setting of the Frank Slide. 


(3) Chapter 6 presents the analysis of the Frank Slide. 


The conclusions of this thesis work are presented in 


Chapter 7. 
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THE STABILEGTY OFOROCKI STOPES 


2.1 INTRODUCTION 


in) the “antroductory chapter the importance of 
establishing the degree of stability of natural rock slopes 
and slopes in open pit mines was discussed. To accomplish 
this, it is necessary to understand which factors influence 
the stability and then assess the stability by developing a 
mathematical model that represents the conditions in the 
Slope. The two types of mathematical models that have been 
developed are the Limit Equilibrium Method and the Stress- 
Strain Analysis. Ultimately, in order to prove whether the 
slope conditions have been correctly modeled and for the 
mathematical model to have any meaning, the field evidence 


must confirm the theoretical predictions. 


The purpose of this chapter is to briefly review the 
analytical methods available for the analysis of slopes and 
to discuss the parameters (necessary for the analysis) 
representing and governing the rock mass behaviour. 
Furthermore, it is the purpose to examine documented records 


of failed rock slopes and to inspect the agreement between 
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the predicted and observed behaviour. 


2.2 ANALYTICAL SLOPE STABILITY ANALYSIS 
Theoretical analyses generally fall into two 
categories, (1) Limit Equilibrium Methods and (2) Stress- 


Strain Analysis. 


Z2e2.1 THE LIMIT EQUILIBRIUM METHOD 


The first step in any Limit Equilibrium Method is to 
assume the shape of a potential slip surface. Having done 
this, Ltt 1s possible; “using *thetequations oigeqiilibriun, te 
determine the shearing resistance required for equilibrium 
to be maintained along the assumed slip surface. The 
required shearing resistance is then compared with the 
available shearing resistance which may have been determined 
from laboratory or field shear tests. This comparison of 
the available and required shearing resistance is. carried 
out by means of the Factor of Safety. The Factor of Safety 
(F.S.) is defined as the value by which the shearing 
resistance parameters must be divided in order to bring the 
potential sliding mass to the point of limiting equilibriun. 
Pavallye ee us procedure is repeated for many possible _ slip 
surfaces, and the most critical one is taken as the one 


controlling the design. 
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The simplest example of a limit equilibrium analysis is 
thea case Oftta Tagid’ block’ sliding* on al= planar surface ”"a’s 
idlustpatedhy inte Faqgure’ 23. 1.8°Ins the absence® of cohesion’ and 
pore water pressure, the activating force down the slope 
would be W ese sina. The force resisting the sliding would 
be Woe coSa@estan-%, where $> is 'the angle “of “frictional 
sliding resistance along the plane. The F.S. would then be 
the activating force divided by the resisting force. In 


equation form 


Past C=ite cosa tettan’ # 
Woe Sina 
= tan 
tan a 


Por l¥aer. SvtGie). 0fttheipoint lof Pimniting *equilibriina, fan- 9 


would be equal to tana. 


invemany soil «slopes; **the “failure” *surface* can be 
approximated by a circular arc. In this case, the potential 
sliding mass is divided into vertical sided slices (Figure 
2.2), and an evaluation is made (for each slice) of the 
terms which express the forces driving the slices down the 
slope and those which resist the motion (Bishop, 1955). At 
the point of limiting equilibrium, the sum of the resisting 
forces and the sum of the driving forces must be equal. 
This type of analysis would also apply to rock slopes where 
theamocks ovproken tup by iclosearregular "joints e*thay at 


may be treated as a soil. 
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From the point of view of satisfying statics, the most 
advanced form of limit equilibrium analysis is the method 
developed by Morgenstern and Price (1965) which is capable 
of dealing with any shape of failure surface (Figure 2.3). 
This ability, of course, is highly desirable and applicable 
to rock slopes where the slip surface is controlled by the 
discontinuities within the rock mass (Piteau, 1970) and 


therefore may not be some simple shape like a circular arc. 


With particular reference to rock slopes, Jennings 
(1970) has developed a method where the failure surface may 
follow several different joint sets (Figure 2.5). With the 
failure surface stepping from one joint set to another it 
may also pass through intact rock. Both the strength of the 
intact rock and the joints are taken into consideration and 


are viewed together as apparent strength parameters. 


In rock slopes, failure may frequently occur in the 
Shape of a wedge (Figure 2.4). The slope and the crest form 
two of the faces with the other being formed by 
discontinuities such as joints or faults intersecting the 
slope. This is a three-dimensional problem and therefore 
must be analyzed as such. This problem is complicated by 
the fact that numerous types of geometrical movements are 
possiblest thativiss toe -Sayperctheviblock-inay)oslidesonetwo 


discontinuities at the same time, or slide on only one _ and 
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break away from the other, and so forth. Goodman and Taylor 
(1967) “treated “this “problem using ‘vectors “while Tonde, et 
al. (1969 LIPO)? and John ™5(1968)-" used stereographic 
projections. Hendron, et al. (1971) have also presented 
analytical and graphical methods to deal with these three- 


dimensional problems. 
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The most logical approach to designing a slope would be 
to determine the state of stress and deformations resulting 
from the stress change associated with the excavation. In 
this way with proper care, the deformations could be 
maintained within tolerable limits and the stresses kept to 
a level below failure. There are, however, several reasons 


why this appreach in general cannot be used (Morgenstern, 


1968). 


In order to perform an analysis where the stresses and 
deformations are calculated requires not only that the 
equations of egquilibrium be known but that the constitutive 
equations governing the material behavior must be 
discernable as well. According to Morgenstern (1968), such 
relationships, except for the probable general form of then, 
are not fully known for rock. The form of these equations 
would most likely be non-linear. The non-linearity together 


with variation in density, pore water pressure in the rock, 
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and the variable boundary conditions presents considerable 
computational “ditivcwlti est Moreover, the analysis would 
only give the change in stress which is not sufficient 
unl ess's the=*inatialiestate Wofl stress? is) known Was well. 
Determination of the initial state of stress is not a simple 


task. 


The computational difficulties have however, to a large 
extent, been overcome in the last decade by the Finite 
Element Method (Zienkiewicz and Cheung, 1967; Desai and 
Abel, 1972). With this method, it is possible to deal with 
variations in density, complex boundary conditions, 
anisotropy and non-linearity with relative ease 
(Guenkvewicz, tet ads, 1970; Stacey;) 19770) exe Nit particular 
reference to rock mechanics, Goodman, et al. (1968a) have 
developed a special one-dimensional finite element to 
Simulate discontinuities within the rock mass. In this way, 
the discontinuities can be assigned different material 
properties than the intact rock, making the analysis more 
realistic. Mahtab (1970) used this type of element in a 


three-dimensional study of jointed rock slopes. 


To the present (1974), however, this numerical tool has 
been used mainly for research purposes. While the Finite 
Flement Method has overcome the computational problems 
inherent in a stress-strain analysis, the material behaviour 


is not generally well understood which renders the analysis 
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CL i imitedavaiuc. For example, for the one-dimensional 
joint-element developed by Goodman, et al. (1968a), three 
parameters are needed to describe the joint behaviour; (1) 
the unit stiffness across the joint, (2). the unit stiffness 
along the joint, and (3) the shear strength along the joint. 
The. .joant stritiness conce pt,.iacconding, .to Goodman ,<!.. «lbs GO 
new 7sohat ~gO »ValUGS.».ATe.-.to, .be. found. in ,weporis and 


publications about joint properties." 


Duncan (1972), in commenting on the use of the finite 
element method for slope stability analysis, expressed the 
opinion that a better understanding of failure and 
postfailure of geotechnical material is needed before 
improvements can be made in stability analysis using the 


finite element method. He goes on to state 


",..while finite element analyses have proven 
useful for studying the occurrence of local 
failure, they have so far been less useful for 


studying slope Sie bi laewy. problems. The 
results achieved so far in using the finite 
method LoL stability calculations have 


provided results which are not significantly 
better than those achievable using accurate 
equilibrium methods of slope stability 
Save vss .cec. 


Undoubtedly, in the future, more research conducted in this 
area will result in a greater understanding of the material 
behaviour making this powerful numerical technique a 


practical design tool. 
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The value of scale models in the study of slope 
stability problems appears to be in the understanding of the 
rock mass behaviour as a whole. For example, Muller, et al. 
(1970) have shown, using scale models, how individual rock 
blocks within a slope can rotate and become unstable. 
Barton (1971a) used scale models to study the effect of 


different orientations of joints with respect to the slope. 


While the models do aid in understanding slope 
behaviour, performing parametric studies using this approach 
would simply be too costly, too time consuming, and _ too 


unrealistic. 


2.4 GOVERNING PARAMETERS IN SLOPE STABILITY ANALYSES 

If .w,examining the«stabilaty of jrockaslopes, tthe tiactors 
which generally must be considered are the geological 
setting, attitude, geometry and spatial distribution of 
discontinuities, ground water conditions and the strength 
parameters of the discontinuities (Piteau, 1970; Patton, et 
Alepbioe WeetOl hese, ufaictons, sf ftabas eeprobabile si thate tite 
greatest uncertainty lies in the values adopted for the 


shear strength parameters. 
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When the shearing resistance is described by the Mohr- 
Coulomb failure criterion, the important shear strength 
parameters are C and d. C is the cohesion intercept and q 
is the angle of shearing resistance. The C value holds an 
unfortunate position in that the F.S. is very sensitive to a 
small change in C (Jaeger, 1971). For example, for a 70° 
Slope, 100 feet high, and $ = 329, the F.S. changes fron 
{PUSS WAThlCe= CPto) 161 i withercioa=text psa (Pentz,2an97d). 
Moreover, this variation in C is well within the range of 


values obtained in shear tests (Jaeger and Rosengren, 1969). 


For many rocks the shearing resistance bears a 
nonlinear relation to the normal pressure (Murrell, 1965; 
Maurer,-© 1\966C2PHObDDS,0 W970)% The non-linear failure 
criterion, however, can be adequately presented as one that 
is piece-wise linear (Morgenstern, 1968). That is to say, a 
particular @ and C are valid only over a limited normal 
stress range. Therefore, even for a non-linear failure 
criterion, the important strength parameters can be 


expressed by a characteristic C and ¢. 


Establishing the § and C parameters from shear tests is 
by no means a simple task. For laboratory testing, the 
problem hegins with obtaining suitable samples. This sample 
problem is discussed further in Chapter 4 in connection with 
the procedures followed in obtaining samples at the Frank 


Slide. 
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Furthermore, there is the question of what size of 
sample constitutes a representative sample. at has 
generally been thought that the larger the sample, the more 
representative at as ‘of. the bfield> Conditions: Formed hus 
reason, Shear boxes have been built as large as 15 x 12 
imchesi*(hoek and #iPentz, 1119.68) Tiend: 115e:8irex? 15.8.) admches 
(Kesmanovici, 2 96y7)c In order to test even larger samples, 
insitu shear tests have been performed in some cases 
(Krsmanovic, et al., 1966; Serafim, et al., 1968; Wallace, 
Sestak, CKIEI)C. While the larger tests may be more 
representative, they are also much more difficult to 
perform. In the case of laboratory testing, larger samples 
are more difficult to obtain and for the insitu tests, the 


cost and time involved often become prohibitive. 


The effect of sample size is further considered in 
Chapter 4 in connection with the shear tests performed on 


the samples from Turtle Mountain. 


A study was made of documented records reporting rock 
slope failures with the purpose of assessing the agreement 
between the theoretical and field behaviour. This study was 
also conducted in order to inspect the range of the strength 


parameters required for a F.S. of unity. A summary of the 
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cases examined is presented in Table 2.1. 


Theas? = avdsiGarrequired for a F.S. of Unity for the 
Failures examined are indeed enlightening. For all the 
cases examined, the 9 is as high as 45° in only one case and 
the lowest $ required is 299°. In the majority of cases, the 
¢@ required is around 300°, The C required is relatively 
small and in many cases is zero. It is interesting to note 
that the § required for a safety factor of unity is very 
Similar to the basic friction angle of 30 + 5 degrees, often 


assumed efor anccha(Gadanyi~siet cahbaz 197 0isicBartonpjalO7 th), . 


For the majority of cases where the shear strength 
parameters have been determined by some type of laboratory 
test, these parameters appear to compare fairly well with 
those obtained from the limit equilibrium analysis. The 
point that must be remembered, however, aS waS pointed out 
earlier, is that small changes in C can have a profound 
effect on the calculated factor of safety. When considering 
the range of cC values obtained from shear tests, the 
correlation is not all that good and in some cases it is 
very poor. For example, in the Twin Buttes slide the 
BAScewithe zero cohesion raforewedge 11 ius /0273. ofliathe dGare 
tinoreased tho 90. 52fepsi acthet »F Sn 1shel..Ceoreandper ft gia pes 
incseased hetorets74 epsipdthethas. fisit. 93h (Seegmiller,;1972). 
Extrapolating this trend to ac of 9 psx, the. lowest value 


obtained from the shear tests, the F.S. is approximatly 7.0. 
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Although this may be an extreme example, it does reveal the 
magnitude by which the F.S. can be in error when calculated 


on the basis of shear tests. 


An item which is of importance and places the majority 
of case histories in a position of doubtful credibility is 
the uncertainty surrounding the pore water pressures at the 
time of failure. If these assumptions are in error, it is 
reflected in the strength parameters required for a F.S. of 
unity; and, in turn, affects the comparison between the 


laboratory and back analysis shear strength parameters. 


2.6 THE RELIABILITY OF ROCK SLOPE STABILITY ANALYSES 


The strength parameters of rock, as previously noted, 
are difficult to determine and, furthermore, the values 
obtained from shear tests may not be the exact values 
operative in the field. The same situation applies to 
establishing the water pressure within the rock mass. The 
PULPOSeeOL the efactoresoitesaretyoMinePan" analysis “1s""to 
compensate for these uncertainties. The numeric value used 
for a F.S. depends on the level of confidence which can be 
placed in the input parameters. This level of confidence 
depends, to a large extent, on the agreement observed 


between the predicted and actual field behaviour. 


During the past several decades, in the area of soil 
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slopes, numerous cases have been documented where a 
relatively good agreement has been observed between the 
mathematical analysis and the field observations. For 
example, for three first-time slides which occurred in 
intact clay, the computed safety factors using the peak 
strength parameters were within 10% of unity (Sevaldson, 
1956; Skempton and Brown, 1961; Kjaernsli and Simons, 1962). 
This type of result has brought about an increased 
confidence in the limit equilibrium method of soil slope 
analysis. With this increased confidence, this method of 
analysis can now be used as a basis for the design of soil 


slopes. 


This type of confidence has not as yet been reached in 
the analysis of rock slopes. In the design of open pits, 
the slope angle is often selected only on the basis of a 
consideration of the local geology and past experience. The 
reason for this is the general lack cf agreement between the 
analyses and observed field behaviour. The F.S. computed on 
the basis of shear test results can be in error by several 
hundred percent. In the Twin Buttes slides, discussed in 
the previous section, the error is as much as 700%. This 
large difference, however, is not always the case. Ate tthe 
Frank Slide (Chapter 6), the F.S. based on the peak strength 
parameters from shear tests on flexural-slip surfaces 
dieters. by only Wugefrom awFisigofeunityslethe jFrankpis lade 


analysis, however, was not as sensitive to a change in 
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cohesion aS is normally the case. Thus the F.S. can differ 
from unity by only a small percentage or by as much as 
several 100%. Since there is a possibility that the 
F.S. based on shear tests may be significantly in error, the 
reliability of analytical rock slope stability analysis has 
not achieved a level where it can be used as the sole basis 
of a slope design. A reasonable amount of judgment is still 


required in the selection of the shear strength parameters. 


Is it possible to improve this position in rock slope 
analysis? Further detailed case studies may improve this 
position somewhat or, alternately, they will further confirm 
the present state of uncertainty. In the case histories 
listed in Table 2.1, as was noted earlier, the strength 
parameters required for F.S. of 1.0 are very similar to the 
basic friction angle as determined on flat rock samples. 
The same situation applies to the Frank Slide, where for a 
P.S.) Of 140, the vaverage ¢@ required (iS 237.29 which is 
identical with the ¢@ measured on a flat sample lapped with a 
eqaese Cigriit (Chapter t6jece Only furthericase Wistories will 
show whether or not this generally applies to all _ rock 
slopes. For this reason it is important to perform careful 
and detailed case studies. Broadly speaking, however, it is 
doubtful if the high level of confidence necessary to design 
slopes to an accuracy of within one degree will ever be 
reached. The material within rock slopes is simply too 


variable and complex to reach this position. 
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Since the reliability of mathematical analyses based on 
insite orjlaboratory Shear tests_.45 gnot? that fhigh, the 
usefulness of these undertakings might be questionable. 
However, the situation exists that for feasibility studies, 
at least, only laboratory-type tests are available on which 
to base the design. Consequently, it is not practical to 
abandon the mathematical analysis unless the design is to be 


based solely on past experience. 


Designing slopes at high factors of safety results in 
low angle slopes which might make an entire open pit 
operation uneconomical. Therefore, it is still necessary to 
work slopes at low factors of safety. Doing this, however, 
necessitates the continuous and careful monitoring of the 
Slopes with the purpose of detecting areas of potential 


instability and remedying them before a catastrophe occurs. 


The situation with natural slopes is somewhat more 
complicated. Since there is no design involved, it isa 
matter of stating whether the rock mass will oor will not 
slide. This matter will be further discussed in connection 


with the Frank Slide. 
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R = the resisting force 


W = weight of the block 


N = normal force on the failure plane 


= W cosa 
The activating force 


The resisting force 


Therefore, the F.S. 


HT 


W sin 
N tan @ 
W cosa 


W cOSsS @ 
W sin 


tan 4 


tan a 


tan @ 


tan 
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Figure 2.1 Limit equilibrium analysissof a sliding block 
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TYPICAL -SLICE 


CIRCULAR 
FAILURE 
SURFACE 


Figure 2.2 [llustration of circular arc failure surface 
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Figure 2.3 Illustration of a general failure surface 
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CHRPLRERe MW! 
THE SHEARING RES ESTANCEH OF 


NATURAL AND ARTIFICIAL POCK SURFACES 
3.1 INTRODUCTION 


Tn rock slope failures, movement generally occurs along 
structural discontinuities within the rock mass. fom, this 
reason, it iS imperative that the factors controlling their 
shear resistance be clearly understood. Numerous shear 
tests have been performed on natural rock discontinuities in 
order to study their shearing resistance behaviour. There 
has generally been no consistency in the reporting of shear 
test data, and the data has been interpreted in a variety of 


ways (Jaeger, 1971). 


This situation has in recent years improved somewhat in 
that the general trend now is to consider the shearing 
resistance along rock discontinuities as arising from two 
components: thee £irst) Mbéingagthestivictional. resistance 
resulting from two flat surfaces sliding over each other and 
the second component arising from the resistance to sliding 


imposed by the geometric irregularities on the rock surface. 
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3.2 THE INFLUENCE OF SURFACE IRREGULARITIES ON THE SHEAR 
STRENGTH 
The frictional component, as determined on flat 


surfaces, can be defined by the expression 


nN 
u 


Nosanspe and.2.2.te.4e8ian.the .ckear. etread 
where S = shear force 

N = normal force, and 

0, = the basic friction angle of the material as 


determined on flat surfaces(Figure 3.1). 


If a sample had a set of regular inclined surfaces 
(teeth) as illustrated in Figure 3.2, and the upper block 
were to slide up the slope as the shearing progressed, the 


sliding resistance can be defined as 
S = W tan (9, + i) ee0ee*ee¢e*eeee#*eesseeeete##eeeeeeee# Se 


where "i" ais the angle of inclination of the teeth upon 
which the sliding takes place. The validity of this 
relationship has been demonstrated by Newland and Allely 


(1957), Ripley and Lee (1961), and Patton (1966). 


The riding of surface irregularities over one another 
can take place only if the sample dilates. When the normal 


stress becomes high enough to prevent dilation, the teeth 
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will be sheared through. When this happens, the shear 
strength will exhibit an apparent cohesion, C, and _ the 


strength envelope can be described by 


S = C + tan b, eeeeoese0evs3#e#0nreeeeeees@eseneteeeeese eee @ ae8 


Combining equations 3.2 and 3.3 to define the shear strength 
over the whole normal stress range results in a bilinear 


failure envelope as shown in Figure 3.2. 


This bilinear law, as demonstrated by Patton (1966), 
describes very well the shear strength of idealized models 
With regular projections. However, for actual rock 
surfaces, there may be a large range of inclination angles 
(Figure 3.3). For this case, the simple bilinear law is 
somewhat unrealistic. Due to the variation in "i" values on 
natural rock surfaces, the change in failure mode from 
Sliding up the slopes to shearing through them will occur 
over a range of normal stresses. Initially, at a very low 
normal stress, there will be only the sliding-up mode of 
failure and eventually at very high normal stresses, there 
will be only the shearing-through mode. In between, there 
could be any degree of intensity of both failure modes 
occurring and, as a result, the failure envelope will be 
curved instead of being simply bilinear. This is supported 
by shear test data which is often best described by a curved 


failure envelope (Murrell, 1965; Maurer, 1966; Hobbs, 1970). 
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Ladanyi and Archambault (1970) realized this and have 


developed a more general shear strength equation than the 


Simple bilinear law. Their equation follows. 
= = : = Co 
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= (oe v ¢f 
where ° = normal stress 
7 = tangential stress 
Vv = the dilation rate due to shearing 
@. = the shear area Tatio; i.e., the fraction of 
the total shear surface, A, passing through 
solid rock 
fF = the average coefficient of friction for the 
contact surfaces of the discontinuities 
7 = the degree of interlocking 
C = the uniaxial compression strength of the 
intact rock 
n = the ratio between the uniaxial compressive 
and tensile strength of rock blocks 


m= (n+1)1/72, 


This equation, while it may conceptually be correct, 
includes parameters like a ,@ ew and’? x which would be 
difficult to measure and for this reason it would be 


difficult to use on a routine basis. 


On the basis of laboratory experiments performed on 
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artificial material with tension discontinuities, Barton 


(1971a) has proposed another law for the peak strength, 
if eS 
o, tan (20 hogy (Caz v5 ) + re) eeecueee 3.4.5 
where C, is the uniaxial compressive strength. 


The simple bilinear law, Ladanyi's equation and 
Barton's equation all included the 6. parameter, and the 
rest of the terms within each equation reflect an attempt at 
including the shearing resistance afforded to the total 


shearing resistance by the surface irregularities. 
3.3 THE EFFECTIVE i ON NATURAL ROCK SURFACES 
Figure 3.3 presents two actual roughness profiles of 


natural rock discontinuities. These traces show the wide 


range in i values which may be present on natural rock 


surfaces. The problem is to determine which is the 
effective i. Patton (1966) was of the opinion that, 
generally, the irregularities could be divided into 
categories of first and second order irregularities. 


Furthermore, Patton and Deere (1971) felt that it is the 
first order irregularities which control the behaviour of 
large slopes and that past histories of tectonic loading and 
weathering will have reduced the interlocking effect of the 


small second order irregularities. 
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Barton (1971b) in examining this question further 
suggests that the $§ +i concept is too simple because it 
assumes that the average angle i remains constant throughout 
the range of normal stresses under which dilation can take 
place. He suggests that the effective value of i depends on 
the magnitude of the normal stress. At very low normal 
stresses, smaller and steeper-sided projections control 
movement. As the normal stresses increase, these smaller 
projections are broken off and the more gentle undulations 
of the first order irregularities control the behaviour. In 
a series of experiments on model material which had been 
fractured in tension, Barton obtained an empirical 


relationship which can be expressed as follows, 
i S=A20 log (CG ) eeeeees¢#*e#e*eeeteensrmeeeseeneoeeee @¢ ol 
we aa 
n 


where C is the uniaxial strength. The assumption used in 
deriving this relationship was that the basic friction angle 


for the material was 309. 


Rengers (1970) and Fecker and Rengers (1971) have shown 
how the effective roughness angle (i) decreases with an 
increase in measuring base length. This suggests that as 
the shear surface area increases, the effective i decreases. 
Therefore, the effective iis a function of the scale of the 


surface area. 
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The approach of determining the 9 in the laboratory 
and then measuring the i component in the field has not been 
applied to practical design of rock slopes. If it has, 
there have been no case histories reported showing that this 
approach can be successfully applied. Therefore, at the 
present time, it is only a concept. This concept, however, 
has greatly aided in providing a rational basis for 
understanding the curvature of failure envelopes which 


usually occur on natural rock surfaces. 


The problem with the approach of measuring the 
roughness in the field is the limited number of suitable 
rock surface exposures with a sufficient areal extent to 
make the measurement. Thisj alse tespecralbly eat rue *satethe 
feasibility stage of an open pit project. Furthermore, in 
order to check this approach and compare it with actual 
failures, it would be necessary to expose the slip surface 
after failure to measure the roughness. To be able to 
examine an exposed failure surface is seldom, and more 
likely never, possible. This limits the applicability of 


this method. 


3.4 BASIC FRICTION ANGLE OF ROCK 


Farlier it was noted that all the recently proposed 


shear strength equations have a parameter included which has 
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heenmaercterred to as the basic friction angle, $,.. This is 
the friction angle as determined from shearing tests on flat 


artificially prepared surfaces (Figure 3.1). 


A fair amount of study has been directed towards the 
problem of how the flat surfaces should be prepared for 


measuring #,. 


Patton (1966) has suggested that the basic friction 
angle be determined by using rock surfaces which are rough 
Sawn. He based his views on the fact that this value was 
Similar to those obtained from field observations of stable, 


unstable and failed slopes. 


courson, (1970,-) 1972) “in following up Patton"s work, 
undertook an extensive investigation to further define this 
basic friction angle for various rock types. Generally he 
found that the frictional resistance depends to a large 
extent on the surface preparation and the surface damage 
which results during the shearing test. The four different 
processes of surface damage which occurred were (1) 
polishing,.)(2); gouging, (3) the generation of rock flour, 
and (4) the formation of an indurated crust. He prepared 
surfaces by using #80 grit, #600 grit and by sand blasting. 
His recommendation was that surfaces prepared with #80 grit 
be used to determine the basic friction angle unless’ the 


field situation warranted a polished or nearly polished 


od 
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surface. 


According to Barton (1971b), the basic friction angle 
can be obtained from surfaces where the shear force 
displacement characteristics show no appreciable peak. The 
granular texture of the rock should be exposed but not to 
the extent that macroscopic interlocking occurs. Sand 


blasting appears to satisfy all these requirements. 


How this basic friction angle, as determined on 
artificially prepared surfaces, compares with the $9 on 
Natural rock surfaces after large shearing displacements, 
still requires some clarification and will be further 


considered in Chapter 4. 


In some of the case histories examined in Section 2.6, 
the $ required for a F.S. of 1.0 was very similar to the ars 
as determined on flat rough-sawn surfaces (Deere, et al., 
1967). This would suggest that the proper way to determine 
b, would be as Patton or Barton suggested. That is to say, 
the flat surface should have a fairly rough texture but no 


large undulations. 


Tt was mentioned earlier that the frictional resistance 
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on flat rock surfaces is dependent on the surface 
preparation. Generally the smoother the surface 
preparation, the lower the frictional resistance (Coulson, 
1970, 1972). A question arising from this is whether or not 
it would be possible to obtain a correlation between the 
microscopic roughness and the frictional resistance. This 
would involve some numerical characterization of the 
surface. Attempts at obtaining this type of correlation 
have been made by Byerlee (1967) on granite and by Myers 
(1962) on metallic surfaces. If such a relationship could 
be established at the microscopic scale, it might also be 
possible to do the same at a larger scale. This idea was 
explored in some detail during the course of this thesis 
work for the limestone rock from Turtle Mountain. The 
remainder of this chapter presents the results of the 


roughness study. 


In order to obtain a mathematical relationship between 
frictional resistance and surface roughness, it is necessary 


to characterize a roughness profile by some numerical value. 


Mechanical engineers in describing the finish on a 
milled surface speak of the root-mean-square value (RMS) of 
that surface. In the British system, a Similar value is 


known as the centre-lineraverage (CLA). inigiact, Snany 
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apparatuses are available which will measure the CLA 


directly and display the result on a dial gauge. 


Numerically the RMS is defined by the following 


equation. (ASAge 1955) = 


N 


RMS = ne y? ax eeeeeee@ee¢eee$ge@eeee#e#e#ee#ee#ee#ee#e#e @ eae 
N 
x = 0 
where N = the number of discrete measurements of the 
amplitude 


y = amplitude of the roughness about a mean of 
Zero 
dx = a constant distance between the amplitude 


readings. 


The CLA value is defined as (ASA, 1955): 


x= .L 
CLA = 1 Ivl dx eeeoeeeeesecsenseeeensve#3ke?eeseenr##erteeeee Bs 
L 
x = 0 
where y = amplitude of the roughness about a mean of 
zero 
L = the distance over which the average is taken 


dx = the length between successive readings of the 


amplitude. 


The difference between these two definitions is in 
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magnitude only. Roughness measuring instruments calibrated 
for RMS will read approximately 11% higher on a given 


surface than those set to read CLA (ASA, 1955). 


Myers (1962) has extended the RMS idea further and has 
given several new parameters for characterizing surface 


roughness. He defined his first new parameter, Z as 


OM 
follows: 


iNebias eeeee7eeeee7eee#eenseoe#eeese#ee @ 3.9 


x = 0 
In words, this is the RMS of the first derivative of the 


profile. 


The second new characteristic, Zee is defined as the 
RMS of the second derivative of the surface roughness. 


Numerically 
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The significance of the Z., parameter, according to 
Myers, is in revealing the roundness of the peaks of the 
roughness profile. Myers correlated these parameters with 


the coefficient of friction on samples of cold-rolled steel 
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disks and found that the best relationship was obtained with 


the Z, parameter. 


Another possible approach to characterizing surface 
roughness, especially on the microscopic scale, is to 
consider the profile to be the result of a stochastic 
process and analyze it on a statistical basis. The 
measurement and analysis procedure for such data has been 
very carefully documented by Bendat and Piersol (1971). The 
majority of information given here with regard to the 
definition of statistical terms was taken from this valuable 


source. 


The four main statistical functions used to describe 
the basic properties of random data are, (1) the mean square 
values, (2) the probability density functions, (3) the 
autocorrelation functions, and (4) the spectral density 
function. The descriptive properties, definitions and 


applications of these functions are given below. 


3.6.1 MEAN SQUARE VALUES 


The mean square value of a sample timerhistory record 
x(t), gives a rudimentary description of the intensity of 
the data. It is simply equal to the average of the squared 
amplitudes of the sample record. In equation form, the mean 


value, MSV, is given by 
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t = N 
MSV = 1 fac dt eeoeoeeoeeoee eee eee ee ee ee oO @ Pe ak 
N 
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Figure 3.4 shows a sample time-history record. The square 
root of the MSV is called the root mean square (RMS). When 
the MSV is obtained with the mean being zero, it is equal to 


the variance. 
3.0. 2 PROBABILITY DENSITY FUNCTION 


The probability density function furnishes information 
about the properties of data in the amplitude domain. It 
gives the probability that the data will assume a certain 
value within a range x and (x + Ax) (Figure 3.4). A special 
case of a probability function is the classical Gaussian 
function characterized by the familiar bell-shaped curve. 
The distribution described by the Gaussian function is 


commonly known as a normal distribution. 
3.6.3 AUTOCORRELATION FUNCTION 


In equation form the autocorrelation function (ACF) is 


described by 
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where 7 is a constant time lag. This function gives an 


indication of the dependence of the values of the data at 
one time on the values at another time. The ACF will have a 
maximum value at zero time displacement. At this point, it 


Will be the same as the variance. 


The principal application of the ACF is to use it asa 
tool to detect deterministic data (e.g. a sine wave) which 
might be masked in a random background. For deterministic 
data the ACF will persist over all time displacements as 


opposed to random data where the ACF will diminish to zero. 


3.6.4 SPECTRAL DENSITY FUNCTIONS 


The power spectral density function of random data 
indicates the manner of the distributicn of the harmonic 
content of the signal over the frequency range from zero to 
Intinity. "if can also be thought of as the amount of its 
mean square value associated with a narrow band of 
frequency, A f. In equation form the spectral density 


function (SDF) is expressed by 
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SDF = [ss (ee (oh erect cmalied ta arty games 
T 
0 


where A, (if) is the Fourier transform of a signal x, (t). 
For a detailed explanation and derivation of this equation, 


the reader is referred to Robson (1963). 


The principal application for the SDF measurement is to 
establish the frequency composition of the data which in 
turn bears an important relationship to the characteristics 
of the system involved. For example, if a system has a 
frequency response function H(f) and a random signal with a 
spectral density function G(f) is applied to the system, the 
output of the system will be a random signal with a power 


spectral density function equal to |{ H(f){2 times G(f). 


When these statistical parameters are viewed together, 
they furnish similar information but in different formats. 
The probability density function furnishes information 
concerning the properties of the data in the amplitude 
domain while the autocorrelation function furnishes 
information in the time domain, and the spectral density 
function gives similar information in the frequency domain 


(Bendatyret ale, 1971) -< 


The physical significance of these parameters will 
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become more evident with the actual digital computation of 


these parameters. 


3.7 SAMPLE PREPARATION 


The first step in exploring this concept was to prepare 
samples with different surface roughnesses and then 
determine the frictional resistance. Withwathisyling mind; 
samples of limestone from Turtle Mountain were prepared for 
direct shear testing by, first of all, using a diamond saw 
to cut specimens from a large rock block. The approximate 
size of these specimens was 1.8 x 1.8 x 3.0 inches. They 
were then cast in Devcon-B epoxy so that the sample would be 
square and fit the shear box exactly. Each half was cast 
separately leaving a space between the halves for the shear 
plane (Figure 3.5). After the casting was completed, the 
two halves were sawn apart using the same diamond saw which 
had been used for the initial trimming. This final diamond- 
saw cut then became the shear surface for the test. Six 


samples were prepared in this manner. 


™n order to obtain different degrees of roughness on 
these shear surfaces, one was left as a diamond-saw cut, 
while another was finished by dry sanding with a #80 grit 
sandpaper, and the final four were lapped with different 
sizes of grit. The sample designation and the corresponding 


surface finish of each surface are presented in Table 3.1. 
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Upon completion of the sample preparation, each sample 
was subjected to direct shear tests using a shear box 
designed and built at the University of Alberta. With this 
apparatus, the normal load was applied with a dead weight- 
lever arm arrangement while the shear load was applied by a 
gear box-chain drive assembly powered by an electric motor 
(Figure 4.5). The horizontal displacement was measured with 
a linear variable displacement transducer (LVDT) and _ the 
Shear load was measured with an 8,00C-pound load cell. The 
shear load and displacements for these particular tests were 
recorded with an X-Y plotter, and the tests were run dry at 
a rate of deformation of 0.05 inches per minute. After 
having completely assembled the apparatus, with the sample 
in place, the horizontal deformation was started and allowed 
to continue until the shear load had more or less levelled 
off. At this point, the normal load was increased and the 
test continued until the shear load had once more levelled 
off. In this way, the sample was step loaded to the highest 
desired normal load. A typical shear load versus 
displacement graph resulting from this procedure is shown in 


Figure 3.6. 


From these load-deformation curves, the shearing 


resistance was noted at each normal load and then both the 
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shear and normal loads were divided by the initial area to 
convert them into stresses. The results for all the samples 


are shown in Figure 3.7. 


The most notable feature concerning the results shown 
IHinkiguren 3.70risr theruwidee variation ogi n friction langbe 
obtained from the various surfaces. From a visual 
inspection and by feeling the surfaces, it seemed that the 
rougher the surface the higher the friction angle. It was 
this observation which initially suggested that it might be 
possible to obtain a relationship between the friction angle 


and the surface roughness. 


The first step in attempting to characterize a surface 
roughness was to obtain a trace of the surface profile on 
the microscopic scale. This was done with a Talysurf 
roughness measuring instrument owned by the Mechanical 
Engineering Department at the University of Alberta which 
produced a graphical output as shown in Figure 3.8. In this 
way, two surface profiles were obtained for each sample, one 


before the shear test and the other after. 


Once the surface profiles were obtained, the next step 
was to digitize them; that is to make discrete measurements 


of the amplitude at regular intervals about some datum. 
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This was done using a digitizer available at the Computing 
Centre, University of Alberta. This apparatus essentially 
converts analog data to digital data and in one of its uses 
converts strip charts to digital representation. The 
digitizer records the X and Y co-ordinates of selected 
points when a cursor is run over the analog signal. The co- 
ordinates are recorded by the digitizer on magnetic tape 
making it possible to handle large quantities of data with 


relative ease. 


For the particular problem at hand (digitizing the 
surface profiles) the apparatus was set to take a reading at 
aMAntervalyot 0.104 anchesi alongil@heaiStrips chart. This 
meant that 100 discrete amplitude measurements were obtained 
for every inch of strip chart. In this manner then, 2,048 
readings were taken for each profile. The reason for this 
particular number is that 2N values are required for the 
calculation of » the “'spectral density function. ™ the 


present context N was chosen to be 11. 


The co-ordinates of each point were initially recorded 
on the magnetic tape in inches with reference to the co- 
ordinate system of the digitizing table. Consequently, all 
the readings had to be converted into micro-inches and also 
translated to some appropriate datum relative to each 
profile. The datum was chosen so that the mean of all the 


Readings oor Giai particular profileM would) be zero. )ithese 
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manipulations were performed on the IBM-360 computer and the 


results stored on a 9-track magnetic tape. 


Later this digital profile data was retrieved and 
replotted using the Calcomp Plotter at the Computing Centre. 
The Calcomp Plotter output for all the surface profiles 
obtained is shown in Figures 3.9 to 3.14. The traces in 
these figures are simply straight lines between the data 
points but the points are close enough to give the 
appearance of smooth curves. These computer plots were 
compared with the original strip charts, and there was no 
visible difference noted between the two. Reproducing these 
profiles from the digital data also aided in ensuring that 


there were no erroneous results in the data files. 


Without the computer, the digitizer and the tape 
storage for information, it would not have been possible to 


deal with such large quantities of data. 


3.10 DIGITAL COMPUTATIONS 

Once the surface profiles were in digital form, it was 
a relatively simple matter to calculate the majority of 
surface characterization parameters discussed in Section 
Beis For example, to calculate the RMS, the requirement is 
to square all the readings, add them, divide by the number 


of readings and then take the square root. This is not 
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mathematically complicated and can be easily done on _ the 
computer even though there were 2,048 readings for each 


profile. 


The only parameter which presented some computational 
difficulties was the spectral density function (SDF). To 
determine the SDF some special techniques are required which 


are worth mentioning. 


There are two methods available for computing SDF's. 
The first is the so-called "Blackman-Tukey" method (Blackman 
and Tukey, 1959) based on computing the autocorrelation 
function (ACF) initially and then calculating the Fourier 
transform tom 8the ACY atThemrourterttransformfoft tthe WCEuMs 
the SDF, and the ACF and SDF are known as Fourier transform 
pairs (Robson, 1963). A newer and more efficient way of 
determining the SDF is the "Cooley-Tukey" method (Cooley and 
Tukey, 1965) which determines the SDF via a fast Fourier 
transform of the original data. Bendat and Piersol (1971) 
give a detailed descripticn of this method and how it is 


used. 


Since the mathematical computations are fairly 
involved, several computer programs have been written to 
calculate the fast Fourier transform of random data. The 
one used during the course of this thesis was obtained from 


the University of Alberta Physics Department and is given in 
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Appendix C. 


3.11 AMPLITUDE DISTRIBUTIONS OF ASPERITIES 

Bowden and Tabor (1956) stated that the results from 
refined profilometer techniques show that the asperities on 
most eee cee have a GausSian distribution. Work done by 
Williamson and Hunt (1968) indicated the same observation. 
It was felt it would be of interest to compare the 
distributions of the surface profiles obtained from the 


limestone samples with the observations of the above- 


mentioned authors. 


From the digital data on each profile, an appropriate 
class interval was chosen and a count was then made of the 
number of readings which fell into each interval. For 
example, the class interval for sample 4 (before shearing) 
was chosen to be 20 microinches, and the number of amplitude 
readings on this sample which had an amplitude range of -120 
to -140 was 65 (Figure 3.18). The results were plotted as 
bars on a frequency-versus-amplitude graph in Figures 3.15 
to P73 720. The theoretical Gaussian distributions are 
Superimposed on each of the bar graphs in order to compare 
the observed with the theoretical. The normal curve 
(Gaussian) can be easily determined by calculating the 


standard deviation of the digital data. 
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Upon examining all those distributions, it appears that 
generally, the actual distribution can be described fairly 
well by a Gaussian distribution. The best agreement was 
obtained for the surfaces lapped with the #200, 400 and 600 
grit before the shear tests were performed (Figures 3.18, 


3219 andesa20)i. 


Probably one of the more interesting features that 
these distributions show is the change which occurs on the 
shear surface during the test. It is most noticeable for 
the samples lapped with #200 and 400 grit (Figures 3.18 and 
3.19). There is a marked increase in the frequency of the 
amplitudes just to the right or positive side of the mean 
while at the same time the highest asperities have 
disappeared. This is quite understandable upon re-examining 
the corresponding surface profiles (Figures 3.12 to 3.13). 
From these profiles it can be seen that the positive peaks 
have disappeared resulting in the flat plateaus. These flat 
plateaus are probably the only areas of contact during the 


friction tests. 


The only distinctive feature about the normal 
distribution curves, as far as characterizing them 
numerically, is the standard deviation which is similar to 


the RMS if the mean is zero. 
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From the digital data, the CLA, RMS, Z, and 23 
parameters were calculated for each surface. The results 
are presented in Table 3.2 together with the friction angle 
obtained from the direct shear tests. The friction angle, 
corresponding with the roughness before the test, was taken 
from the initial slope of the shear versus normal stress 
curves (Figure 3.7). The results tabulated in Table 3.2 are 
plotteds ine Figurec3.21itov3i23i “Since thesCLArand RMS sare 
Similar, except for the 11% difference in humerical 
magnitude, only the CLA values were plotted against ¢. The 
best correlation is between # and the Zo parameter. The 
coefficient of correlation between § and Z, is 0.881. Ina 
statistical sense, the required correlation coefficient, at 
a Significance level of 1%, should be greater than 0.708 for 
there to be a linear correlation. Since the correlation 
coefficient is higher than this, it can be stated that there 
is a linear correlation between Z, and @. The same is true 
for Z, but not for the CLA parameter. The correlation 
coefficient between ¢ and the CLA parameter is only 0.674. 
This type of result is very similar to the results obtained 
by Myers (1962). He obtained a correlation coefficient with 


Z, of 0.84 and only 0.54 with RMS. 


A typical autocorrelation function (ACF) plot is 
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presented in Figure 3.24. Of all the ACF plots, there was 
really nothing distinctive about them which could be used to 
depict a single characterizing parameter. Therefore, this 
function could not be used to define a parameter with which 
the § could be correlated. The value of ACF lies in the 
fact that it is a mathematical tool which can be used to 


determine a harmonic wave masked in the random data. 


As far as establishing a single roughness parameter, 
the same comments apply to the spectral density functions as 
they do to the ACF. A typical SDF is presented in Figure 
3.25. The SDF does not furnish much more information than 
the ACF except that the SDF is presented in the frequency 


domain. 


Since the spectral density function has been used to 
shgay roughness profiles of airport runways and highways 
(Houbolt, 1961; Hutchinson, 1965), it is worthwhile to see 
how it has been applied in the transportation field. Figure 
3.26 gives the criteria as presented by Houbolt (1961). It 
should be noted here that the ultimate purpose of spectral 
density measurements in roadways and airport runways is 
somewhat different than simply determining a Single 
characterizing value for a particular roadway roughness. 
Their interest is more in determining which wavelengths (or 
frequencies) will cause excessive vertical motian to be 


induced when a vehicle passes over the road (Quinn, et al., 
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V970)% The roughness spectrum indicates which wavelengths 


contribute to the total pavement roughness. 


The correlation of roughness parameters with the 
friction angle generally follows a linear trend with @ 
increasing as the roughness parameter increases. An 
increase in the roughness parameter indicates an increase in 


surface roughness. 


The main reason for the amount of scatter obtained in 
the correlation data is that various types and degrees of 
surface damage occur during the shear test. The type of 
damage greatly affects the shearing resistance. Two samples 
with the same surface finish may exhibit a different 
shearing resistance because on one there may be a lot of 
gouging and on the other, there may be surface polishing. 


This is especially true for the smoother surfaces. 


The correlations of roughness parameters with ¢ are 
good enough, especially with 2Z,, that upon knowing the 
roughness parameters, a fairly good estimate could be made 
of ¢. This approach, however, is not practical since the 
time and effort involved in calculating the roughness 


parameter is far greater than Simply doing a shear test. 
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The advantage of performing a study like this at the 
microscopic scale is that it is possible to obtain a 
roughness profile which is representative of the type of 
surface. That is to say, a profile like the one shown in 
Figure 3.10 could have been taken anywhere on the sample and 
it would not have been significantly different. With the 
representative profile, it is possible to develop the tools 
necessary to define the roughness and establish what 
parameter is most clcsely correlated with q. Once this has 
been done, there is conceptually no reason why the same 
principal could not be transferred conte field scale where 
it would be simpler to determine the roughness parameter 


than to perform a shear test. 


The problem with testing the developed nathematical 
tools in the field is two-fold. Firstly, it would be 
necessary to have rock discontinuity exposures in the field 
on which a representative field roughness profile could he 
obtained. Three surface profiles over a five-foot distance 
(Figure A.7) are shown in Appendix A. Certainly these could 
not be considered as representative because the larger 
wavelengths are not complete. This means that the roughness 
profile would have to be obtained over a much larger 
distance; i.e., 100 feet. These types of exposures seldom, 
if ever, exist. Secondly, the associated shear tests at 


this scale could never be performed. 
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TABLE 3.1 


SAMPLE DESIGNATION FOR ROUGHNESS STUDY 


Sample Number Surface Finish 
Af Dry sanded 
2 Diamond-saw cut 
3 Lapped with 45/80 grit 
4 Lapped with #200 grit 
e Lapped with #400 grit 


6 Lapped with #600 grit 
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SHEAR STRENGTH,S 


Figuirews oL 


NORMAL LOAD, N 


Failure envelope for specimens with flat 
surfaces 
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SHEAR STRENGTH, S 
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Figure 3.2 Failure envelope for multiple inclined 
surfaces 


NORMAL LOAD, N 


a) Second -Order 
Irregularities 


i=13° b) First-Order 
Irregularities 


Approximate Scale: ! inch = I foot 


Figure 3.3 An example of a natural discontinuity 
iilustrating first— and second-order 
irregularities (after Patton, 1966) 
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Figure 3.4 An example of a time-history record x(t) 
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LIMESTONE SAMPLE 


STEP | 
MOLD 
DEVCON- B EPOXY 
SET DEVCON- B 
SHOEI ya LIMESTONE SAMPLE 
MOLD 
DEVCON-B' EPOXY 
DIAMOND-SAW CUT 
STEP 3 


Figure 3.55 Casting of “limestone samples 


a -noovaa Tae 


guamae 3M0T@AMII— 
aJjom | 
YxONS @-~noovad 


s 


TUD «waAe-aMo 


61 


G2:0 


qOTd uot}ewtOFEep snstaA peotT Aeaeys ‘papeoTt days TeorddAA wv 9°¢€ eanbty 


(SSHONI ) NOILVWYO0IS3QG 
OS:-O0 Gc-O 


ip) 

ae 

Mm 

Pp 

OSz » 

% 

= 

S87 NI QvOl TvWNWYON——~ 062 oO 
0901 00S 
z 

(98) 

are) 

OS Z 


NI 1-€ :V3uNV 


3 ; 
2261 ‘> 934 


b ‘ON JATIdWVS 


s 


=i 


o80i 


gata JAMAOW——+ cet _— | a : z 
| | 1s 


Bie 
Saé ; i : ao = =. - a 


ad 
= a So 
> 


oe (23HOvi)} vOITAMAOARS 


& toiq motssnicish aveisv bsol usefe .bebsoi qete Isciqy? A a.€ sauprs 7 = 


62 


sessouyhbnor 
SNOTIPA YIM SOoeFIANS YOoT JoZF sadopToeause seiantteqg L°¢ sanhty 


(1Sd) SS3YLS T1VWHON 
009 OO 002 


Ww) 
a 
OG Olesae 
> 
A 
eyere ac 
ys) 
Mm 
WY) 
Ww) 
LIU9 O009¢0 Oo¢ 
LIYD OOb#y = 
LIND O02# o 2 


LIYD O8/GSbty# a 
IND MVS-QNOWVIGs O00Ov 


Q30NVS Aude 


TIND OB\eS* ¥ 


Es. | Tina O0S% © 
Tind OO ** : 
Tind OOa* 4 : 
= oe 
008 oF | 005 . oe 
os (729) e2ea3ATe JAMAON ae ae 
_ avolisy dtiw aesstiua Avo1 102 esqofovas saviiel ‘-& espera 
‘B3829nnpuo2 _ 


63 


ie a ee ee ee oe ees ae Fa eee ae Pee ede. eae eee 


MADE IN ENGLAND 


Ryn eeepc preseye epee eepe ee ge a ese = a 
BY ee at at ee era ng coe ae race eaten eed Ay Vee eel sy goog deer ee ey ee ree 
2 ee ee A aie . iS sok es : - = ee 


Paeroe ae = = “hse = eee 


1 URN CESS een UL cone ee Set RO Soin: AIS MRE Min core Br: 
Ho ee te 


Se cay Re page PN ete wy) ; Pn hes u ar ea 
po i ae eres Ra rar ane SAE pnt Sree tN gle ny eset Pain Wee hee = ls oe ee ee Ne 
rf oe ao ‘ ; ‘ st SS Sie 
= ae ee SEAS ee ans at Spee ean Suter aes me ess pee ehas Peg ot USE eae. Be eae nanticmann = Btosese 5S, bra 9 a 


Ape ae Sees ! 1 H 
SS eM ero ae eer i i eee File esp meena Fete ens et) pee eee el Se Sd eg aes 
upd ; \f is 


Ge rT] Fo SN = eS ene na tc Pn eS os gs 
Pa) ee 5 aE Pe pee earner ay SS Ae Cee SS ee ee ee = os 
Aetetoe feat = ee — Bia een = 


a a ee oO 


{ 


ee oe ee ee ice ee ae, et I Me 


Pe 
he 
| 
i 
j 
t 
4 
{ 
y 
tj ei] 
} 
ee 
| 
i 

i 

| 

| 


SCALE: VERTICAL [| DIve= 10G@ 4” 
HORIZONTAL | DIV= 0-002" 


Figure 3.8 An example of a surface trace from the 
Talysurf 
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Figure 3.15 
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Amplitude distributions of asperities for 
dry sanded sample 
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Figure 3.16 Amplitude distributions of asperities for 
diamond-saw cut surface 
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Figure 3.17 Amplitude distributions of asperities for 
surface lapped with #45/80 grit 
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Figure 3.18 Amplitude distribution of asperities for 
surface lapped with #200 grit 
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Figure 3.19 Amplitude distribution of asperities for 
surface lapped with #400 grit 
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Figure 3.20 Amplitude distribution of asperities for 
surface lapped with #600 Grout 
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Figure 3.26 Limiting spectral criteria for roadways 
(after Houbolt, 1961) 
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CHAPTER IV 
SHEAR TESTS ON NATURAL ROCK 


DISCONTINUITIES FROM TURTLE MOUNTAIN 


4.1 INTRODUCTION 

In back analyzing a slide the factor of Safety is taken 
to be unity and in doing so, it is possible to determine the 
Shear strength parameters required for this condition. The 
main purpose of the peek analysis, however, is to determine 
how reliable the analysis is and to see af the safety factor 
of unity could have been predicted by determining the shear 
strength parameters from laboratory shear tests. Thus in. 
the context of the Frank Slide, it was vital that the shear 
strength along the failure surface be determined from 
laboratory shear tests on samples which contained 
discontinuities representative of the failure surface. This 
Chapter presents the methods used to obtain the necessary 
Samples, how they were tested and the results obtained on 
the various types of discontinuities. Included as well is a 


discussion on the implications of the test results. 
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4.2 SAMPLING LOCATION 

Upon exploring Turtle Mountain and the surrounding area 
in some detail, it was Soon realized that due to the 
magnitude oe the slide and the rough terrain, the only 
accessible place to obtain Samples would be in the rock 
debris atong bthekeroade A secondary road to the west of 
Highway 3 Passing through the rock debris Was ~found® te” “he 
the most suitable location. The major portion of the Slap 
surface at the Frank Slide was in the Livingstone Formation 
aS was the rock mass which “SLid. Therefore, the rock 
anywhere within the slide debris would be representative of 
EhespecekAvinvolvede ind the Slide. A description of the 


Livingstone Formation is given in Chapter 5. 
4.3 SAMPLING PROCEDURE 


After considering the geology of “Turtle Mountain "iH 
was decided to make an attempt ae saul iii three different 
types of discontinuities; namely, joints, bedding planes and 
surfaces which showed Signs of previous shear movement along 
them. The latter type will be referred" to ‘ash!*"flexural- 


Slip" surfaces. 


The boulder chosen ‘as. the Sample block is shown in 
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Figure 4,1. The bedding planes and joints can be very 
clearly seen in this Figure. Figure 4.2 shows a flexural- 


Slip surface within the sample block. 


Two methods were used to obtain samples from =this 
boulder. One was to simply sever smaller blocks from the 
boulder by wedging a chisel into joints and bedding planes 
which were slightly open. This method worked fairly well 
for obtaining blocks which contained joints. However, since 
the bedding planes were spaced approximately two feet apart, 
it was impossible to obtain a single block with a bedding 
plane passing through it. The block would have Simply been 
too large. Therefore, to use this method of obtaining 
blocks which contained bedding planes meant that each side 
of the bedding plane had to be worked at separately. Some 
bedding plane samples were obtained in this way; however, 
the number obtained suitable for testing was low in 


comparison to the number of attempts. 


Generally, the problem with using only hand tools and 
prying away the blocks was that the size of the blocks was 
too large to handle so that when they were severed from the 
large boulder, they would fall to the ground and shatter. 
From the shattered rock it was possible at times to obtain 


Small samples suitable for preparing a lab specimen. 


The other method used was to core around desirable 
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discontinuities with a masonry drill manufactured by Delro 
Industries Ltd., Winnipeg, Manitoba. This machine is 
IRtend claw fores drilling large diameter holes in masonry, 
reinforced concrete or other hard material. The diameter of 
the diamond core barrel was ten inches. Figuye 4.3 5 shows 
the coring machine anchored to the sample block at Turtle 
medneaa a while Figure 4.4 shows a core obtained in this 


Mann ear. 


THEY ecoring “method) sic advantageous from the point of 
view that the total volume of rock to be handled is more 
manageable and the surface area of the discontinuity is 
fairly large. This procedure, however, does have some 
severe problems associated with it which limits its uses. 
The biggest problem is the retrieval of the core. 
Sometimes, it was fortunate that a discontinuity existed at 
the maximum depth of coring (approximately 12 inches) and 
then there was no. problem. However, if this was not the 
case, the core recovery was extremely difficult and 


sometimes impossible. 


In most cases, whether it was a core or block sample 
being worked on, the specimen would fall apart into several 
pieces. In this case each piece was marked and the block 
reassembled and then wired together for transportation back 


to the lab. 
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Besides the Samples which contained natural 
discontinuities, solid blocks of limestone were brought back 
to the laboratory from which the samples having artificial 


surfaces would be prepared. 


From the larger blocks and core samples, smaller 
Specimens were cut to approximate size using a diamond saw. 
With this saw, it was impossible to get the sample exactly 
Square and in perfect shape for the shear box. Tor “this 
reason the samples were brought up to final size by casting 
them in Devcon-B plastic steel epoxy. It was done in a 
Similar manner to that mentioned in Section 3.7. This 
casting ensured a snug fit in the shear box and also 
peuaeed in a flat horizontal surface for seating and 


loading the sample. 
4.5 TESTING APPARATUS 


Upon completing the sample preparation, each of the 
samples was subjected to a direct shear test. Two types of 
Shear boxes were used for this purpose. The first, designed 
for 2" x 2" samples, was built at the University of Alberta 
and.) is Wshown@iin®s igure? 825. The shear load with this 
apparatus is applied by a gear box chain drive assembly 


powered by an electric motor. The normal load is applied 
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with a dead weight lever arm arrangement. Te horizontal 
and vertical deformations were measured with linear voltage 
displacement transducers (LVDT) while the shear load was 


measured with an 8,000-pound load cell. 


The other machine used was a ten-ton Wykeham-Farrance 
direct shear machine shown in Figure 4.6. The shear load 
with this machine is applied through a variable gear 
arrangement powered by an electric motor. The normal load 
is applied through a loading yoke by a hydraulic ram. [fn 
order to maintain a more accurate control over the normal 
load than by simply using a dial gauge on the hydraulic 
pressure, a load cell was incorporated between the loading 
yoke and the load cap. Strain gauges on the load cell were 
connected to a Budd strain indicator which was then used to 


measure the normal load. 


With the Wykeham-Farrance machine, the shear load was 
measured by strain gauges which had been mounted on the ten- 
tong provang) ang. AS with the small shear box, the 
horizontal and vertical deformations were measured with 
LVDT's. To measure the vertical deformation, an LVDT was 
mounted on each corner of the top plate. With these four 
LVDT's, it was possible to measure the rotation of the top 
half of the sample and the top loading plate. The apparatus 
was originally designed for 12" x 12" samples. But this was 


modified with spacer blocks so that 10", 8" and 6" square 
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Samples could be tested as well. 


With both shear boxes, the shear load and horizontal 
deformation were recorded with an X-Y plotter. The vertical 
deformations were recorded with a Hewlett-Packard digital 


recorder. 


In order to obtain surface profiles of the natural 
discontinuities, a small frame was built on which an LVDT 
and stylus arrangement was mounted. The vertical movement 
of the stylus was sensed by the LVDT and the output from the 
LVDT was recorded directly on an X-Y plotter. Figure 4.8 
Shows this apparatus. The X-Y¥ plots were later digitized 


and then replotted with the Calcomp Plotter. 


The shear testing program was performed in two stages. 
To begin with, five series of 6" x 6" Samples (nominal area) 
were prepared and tested while later, in the second stage, 


the 2" x 2" specimens were tested. 


In the first five series of samples, three series 
contained natural discontinuities while the other two had 
artificially prepared surfaces. An example of an artificial 
surface would be a diamond-saw cut. The natural surface 


Studies were composed of a series with bedding planes, a 
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series with joints and a series with flexural-slip surfaces. 
The flexural-slip surfaces were prepared so that the 
shearing would be parallel to the striations. Figure 4.7 
shows typical samples of these series. Of "ther rtwo 
artificial series, one was lapped with a #45780 grit and the 
other was made with a diamond saw. Each series contained 
four samples except in the case of the bedding plane series 


which contained only three samples. 


The first series of samples was then tested in the ten- 
ton Wykeham-Farrance direct shear machine. The tests. were 


run dry at a rate of deformation of 0.048 inches/minute. 


The general procedure after the machine and sample had 
been assembled was to apply the normal load and then 
initialize the recording devices. When the set-up was 
completed, the horizontal deformation was started and 
continued until the total deformation was approximately 0.8 
inches. At this point, the normal load was held constant 
and the direction of the shear displacement reversed and 


sheared back to its initial position. 


Upon removing the samples from the shear box after the 
test, it was noticed that the actual area of contact had 
been very small. Figure 4.9 shows this remarkably well. 
Since the number of samples obtained in the field was 


limited, it was decided to make further use of these 6" x 6" 
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Samples by cutting 2" x 2" samples out of them where there 


had been no contact. 


Thus, in the second stage, two sets of 2" x 2" samples 
were prepared for testing in the small shear box. These 
small samples were prepared only for the bedding planes and 
flexural-slip surfaces. Each of these two series contained 
four samples. After the preparation was completed, they 


were tested in the same way as the larger samples. 


For each of the shear tests run, a shear load versus 
deformation plot was obtained directly on the X-Y plotter. 
These X-Y plots were later digitized and replotted on the 
Calcomp Plotter. This was done so the curves could be 
plotted on a more appropriate scale and in § a more 
presentable form. A typical result of this type of plet is 


shown in Figure 4.10. 


The digital vertical deformation data was punched up on 
computer cards and then reduced on the computer and plotted 
on the Calcomp Plotter. Figure 4.11 shows a typical 


vertical deformation curve. 


From each of the load-deformation curves, note was made 


of the peak strength and the post-peak strength in both 
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directions. The post-peak value for shearing away from the 
initial position was chosen at a point after the shear 
resistance had essentially levelled off shown at Point "a" 
in Figure 4.10. The reverse post-peak value was then taken 
at a point when the sample had returned to the same position 
as it was at Point A (Figure 4.10). Since there was a fair 
amount of slack in the machine, this meant that the reverse 
post-peak value was read at Point B, to the left of Point BR. 


This is shown in Figure 4.10. 


The slack in the shear machine also showed up on the 
vertical deformation curves. In Figure 4.11, if the reverse 
portion of the curve were displaced to the Sightyipyesrhe 
amount of the slack, it would fall much closer to the 
vertical deformation curve when Shearing from the initial 


position. 


The shear loads and corresponding normal loads were 
then divided by the initial area to convert them to stresses 
and then plotted as shown in Figures 4.12 to 4.16. The 
failure envelopes drawn through the data points on these 
figures were obtained from a linear regression analysis. 
The only case where the linear regression analysis was not 
done was on the diamond-saw cut surface. Table 4.1 gives a 
summary of the strength parameters obtained on each type of 


discontinuity. 
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AS was mentioned at the beginning of the chapter, the 
main purpose for performing the shear tests on the natural 
discontinuities was to determine strength parameters for the 
Prank Slide stability analysis. There are, however, some 
interesting observations which have emerged from the shear 
tests which are worth mentioning in connection with the 


interpretation of the data. 


4.8.1 THE ULTIMATE FRICTION ANGLE 


In Chapter 3 a discussion was given on the Significance 
Of thepbaste frictiontangle, Oh» and the question was raised 
/ as to how it compares with the friction angle on natural 
rock surfaces after large shearing displacements have 
occurred. This friction angle on natural surfaces after 
large shearing displacements will be referred to as the 


ultimate friction angle, ¢@. . Determining the @, on natural 


U u 


rock surfaces is complicated by the fact that natural 
surfaces are seldom, if ever, flat thereby making it 
difficult to isolate the resistance due to $, and that due 
to the geometric component. For this reason, it would be 
desirable to measure fy on flat artifically prepared 
Surfaces if it could be shown that the Qieon Lhatucal 


surfaces is the same as §,. 
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In order to measure q on natural rock surfaces, itis 
necessary to isolate the §, term from the geometric 
component of the shearing resistance. This can be done by 
considering the post-peak shearing resistance, when shearing 
away from the initial position, as being d,+ i and when 
Shearing in the reverse direction, back to the initial 
positron, was bur i. The average of the post-peak shearing 
resistances in both directions will be $,. To eliminate 
the i component in this manner requires that the i value is 
the same in both directions. This was considered, as 
discussed in Section 4.7, in selecting the post-peak values 
from the load-deformation curves. Furthermore, this also 
implies that shearing in the reverse direction should only 
be continued until the sample reaches its initial position. 
To shear past the initial position would be like Starting a 
new test since the shearing surface would then be on the 
Opposite side of the irregularities. Tuy yhact,) ebyaadeing 


this, it is possible to do two tests on one Sample. 


The ultimate strength parameters quoted in Table 4.1 
and the ultimate strength envelopes drawn in Figures 4.12 to 
4.16 are based on the averages of the post-peak values in 


each direction. 


It should be noted that even on flat surfaces like the 
lapped surfaces, there can be a slight difference in 


Shearing resistance in each direction (Figure 4.15). This 
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points out that there are other factors which can cause 
minor directional shearing differences. However, if there 
is a substantial difference in each direction, the major 


contributing factor would be the geometric component. 


Taking an overall view of all the ?, values obtained 
(Table 4.1), it is noted that the $b, LOeeSrockOhsurtaces Vis 
not aunigue value. For the tests performed, the a varies 
anywhere from 37.29 on the lapped surfaces to 14.0° on the 
flexural-slip surfaces: The main factor which seems to 
control the #, value is the type of surface damage which 
occurs during the tests which, in turn, depends on the 
condition of the shear surface before the shear test. On 
the flexural-slip ‘surfaces ‘and ‘joints, a shiny indurated 
crust had formed at the points of contact which was hardly 
noticeable on the bedding planes. On the bedding planes 
there was more gouging than on the other two natural 
surfaces. The shiny indurated crust on the joints and 
flexural-slip surfaces are most likely responsible for the 


low ultimate strength. 


The shear test results show that the ¢, as determined 
by a reversal shear test on natural rock surfaces can differ 
Significantly from the 6, asiideternined. ont flat @artitical ly 
prepared surfaces. This is contrary to the assumption often 
made which states that they are the same. This is at least 


true for the limestone within the normal stress range at 
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which the shear tests were performed in this study. Ata 
very high normal stress, it is conceivable that the 
indurated crust would break and then the ultimate shearing 
resistance would be the same on all the surfaces. Although 
this may be the case, the shear behaviour at the high normal 
stresses is not applicable to the understanding of rock 
Slope failures since the normal stresses along the _ slip 


surface in rock slopes are generally relatively low. 


The shear strength of clays after large shearing 
displacements is defined as the residual strength. It is a 
unique value dependent only on the type of material and to 
some extent on the normal stress. The shear strength of 
rock surfaces after large shearing displacements is 
dependent on the same variables but it is also greatly 
affected by the surface roughness. Therefore within context 
of the the definition of the residual strength for clays, 
there is no residual strength for rock surfaces. FO Ps pthais 
reason the shear strength on rock surfaces after large 
shearing displacements will be referred to as the ultimate 


strength. 


Whether the $,, as determined on natural surfaces from 
a reversal test, has any significance in the field will only 


become apparent from case history studies. 
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428e2 THE EFFECT OF SIZE .OF SAMPLE ON THE SHEAR STRENGTH 


The shear tests performed in thas investigation 
included samples which contained the same type .of 
discontinuity but were different in size. The larger 


Samples had a nominal area of 6" x 6" and the smaller ones 
were 2" x 2", It was explained earlier how the 2" x 2m 
Samples were obtained from the larger samples, This 
procedure had the added advantage of being able to 
investigate the size effect on the shearing resistance. 
When examining the results in FAGURGSEH 1 2eandsd. G3) bie iS 
hevedesthatpetheressis no Significant difference between the 


two sizes outside of the usual scatter. 


This result is probably not that SUP PEISING wet ean 
examination is made of the surface profiles of these 
Samples, Figure 4.17 shows three typical profiles for the 
bedding planes and Figure 4.18 shows the profiles for the 
flexural-slip surfaces. These profiles were intentionally 
plotted at the natural scale for the most realistic 
comparison. From Figures 4.17 and 4.18, it can be seen that 
the scale of roughness on 2" x 2" Samples would not be 


Significantly different from that on 6" x 6" samples. 


A similar result was observed at Imperial College where 
tests were performed on a 100-ton shear machine with a 


Sample sizeof. 15" x. 12" and as .well on. a small ten-ton 
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field shear machine, The two machines gave Very Similar 


results (Hoek, NOT Va'ye 


The shear tests performed in this thesis and the ones 
reported in the literature indicate that the shear strength 
is not influenced by the size of the Sample as long as the 
Scale of roughness remains constant. if the Scale or 
roughness differs on each Sample size, as shown in Figure 
4.19, then the shear strength would be different on each 
Sample because of the different geometric effects, At the 
scale of ‘hundreds of feet, it would never be possible to 
perform a shear test, Therefore, at this large’ scale, it is 
hecessary to determine the Shearing resistance at a much 
smaller scale and then adda to Pt" "the effective 1 as 
determined by field roughness measurements. The” “inportant 
point to remember is that the shear strength is cnly 
affected by the size of the Sample if the scale of roughness 


changes with the size. 


4.8.3 THE SHEAR STRENGTH OF BEDDING PLANES AND FLEXURAL- 


SLLP SURFACES 


The shear strength envelopes for the bedding planes 
(Maguire Yet25n care considerably different than those for the 
flexural-slip surfaces. Tivs' "at firet glance, vo 9 somewhe < 
of an anomaly since the flexural-slip surfaces are the 


result of shear movements along bedding planes. Therefore 
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it would seem logical that the flexural-slip surfaces would 
not exhibit a peak strength and the ultimate strength would 


be the same as on bedding. 


This, however, is not the case. The reason EOe Chia Ss es 
that the mode of shear failure in the field was much 
different than that which is observed in the laboratory. 
The shearing displacements which occurred along the bedding 
planes during the folding cf the strata occurred at very 
high stresses and temperatures. These conditions, together 
with the geologic time scale of shear movement, allowed 
recrystallization to take place. The recrystallization 
resulted in the shearing surfaces remaining in intimate 
Gontact. By subjecting a bedding plane to a shear test in 
the lab, the mode of shear failure and the resulting sheared 
Surface will be much different than that which occurred 
during the flexural folding. Thus, although the flexural- 
Slip surfaces are the result of flexural slipping on = the 
bedding planes, their shear strength behaviour bears no 
direct relationship to the bedding planes in the realm of 


laboratory shear testing. 


In the previous section it was pointed out that the 
shearing strength of natural rock surfaces is a function of 
the initial surface roughness and the type of surface damage 
which occurs during the shear tests, Bothy the ana tiad 


surface roughness and the shearing surface damage differ on 
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the bedding planes and the flexural-slip surfaces; 
therefore, it is understandable that the shear strength will 


be different. 


TABLE 4.1 


SUMMARY OF SHEAR STRENGTH PARAMETERS 
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Peak Ultimate 
Type of Sample @ C* (psi) @ CG (psi) 
Bedding plane Die! 38 32.3 8 
Flexural-slip 
surface 28.20 ae 15.6 18 
Joints 52.0 25 14.0 deo 
Diamond-saw cut 200 0) 297.0 ) 


Surface lapped 
with 45/80 grit ---- -- Sire 5 
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Figure 4.1 The sample block within the slide 
debris 


Figure 4.2 A flexural-slip surface within the 
sample block 
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Figure 423 (€oringjon)the samplesblock 


Figure 4.4 A ten-inch diameter cored sample 
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Figure 4.5 The modified shear box for the two- 
inch samples 


Figure 4.6 The Wykeham-—Farrance shear box 
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Figure 4.9 A joint surface after shearing 
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Figure 4.12 Failure envelopes for the bedding planes 
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Figure 4.13 Failure envelopes for the flexural-slip 
surfaces 
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Figure 4.14 Failure envelopes for the joint surfaces 
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Figure 4.16 Failure envelope for the diamond-saw cut 


surfaces 
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Figure 4.17 Natural surface profiles of a bedding plane 
(sample B-3) 
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Figure 4.18 Natural surface profiles of a flexural-slip 
surface (sample S-2) 
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CHAPTER V 


THE GEOLOGY OF TURTLE MOUNTAIN 


5.1 INTRODUCTION 

In 1903 on April 29, at 4:10 AM, a gigantic rock slide 
occurred on the east face of Turtle Mountain. The slide 
destroyed the southern end of the town of Frank, a small 
coal mining town in the Crowsnest Pass area of southwestern 
Alberta (Figure 5.1). The Canadian American Coal Company, 
for whom most of the inhabitants of Frank worked, was at the 


time mining a coal seam at the base of Turtle Mountain. 


The slide debris moved down the east face of the 
mountain across the entrance to the Frank mine, the southern 
end of the Town of Frank, the main road from the east and 
the Canadian Pacific main line through the Crowsnest Pass. 
The slide then continued up the opposite side of the valley 
before coming to rest approximately 400 feet above the 


valley floor and about a mile from the base of the mountain. 


Approximately 70 people lost their lives in the slide 
and several miners were trapped in the mine for a 


considerable time. The disastrous results of the slide 
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immediately aroused the question in the survivors! minds, 
"will another slide occur and destroy the remainder of the 
town?" This fear resulted in several investigations to 
determine the conditions of fTurtle Mountain and the 


possibilities of another slide. 


These investigations, considered in detail below, 
presented conflicting reports as to the geological structure 
of the mountain. In order to correctly establish the 
geology, several trips were made to the slide and an 
examination made of the area. This chapter gives the 
geology of Turtle Mountain as envisioned from the 
observations made on the visits to the apnea wand’ von Fhe 


existing documented information. 


5.2 HISTORY OF INVESTIGATION 

Immediately after the slide in May, 1903, an inspection 
was made by R. G. McConnell and R. W. Brock of the 
Geological Survey of Canada. Their report, completed within 
a month of the disaster, gave a general survey of the 
geology of the mountain. They concluded that the slide 
occurred across, rather than along, bedding planes ana 
believed that the primary cause for the slide was to be 
found in the structure of the mountain. te Wase thelr 


Opinion that if there was any further danger of another 
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Slide, it would come from the north peak. 


Brock visited Turtle Mountain again in 1909 and 1910 
and became increasingly concerned about the possibility of a 
Slide from the north peak (Brock, 1910, 1911). POWs al result 
of his report in 1910, a commission was appointed to study 


the conditions of Turtle Mcuntain. 


The’ report of this commission (Daly, ie ae ele) 
again emphasized the danger of another rock slide from the 
north peak of Turtle Mountain. Their reasoning was that the 
structure of the mountain is the same along all profiles 
from the south side of the south peak to a point well past 
the north peak. They concluded that because the north peak 
is similar in structure to the portion which fell away, 
Sliding was likely to occur again at the north peak. It was 
in this report that a cross-section (Figure 5.2) showing the 
geology of the north peak was published. This cross-section 
has appeared in many publications, notably texts by Coates 
(1967), Leggett (1962), Leggett (1973) ,°o9 Longwellyi Pet© aie 
(14269) ,AlReynoldss) (1961) Ft eSharpesm(sssjsenzarubay Veal: 
(1969), and in reviews by Terzaghi (1950) and Ter-Stepanian 


(1966). 


The structure of Turtle Mountain was described as a 
monocline of Palaeozoic limestone dipping to the west about 


fifty degrees. The limestone which formed the northerly 
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trending ridge of the Blairmore Range, of which Turtle 
Mountain is a part, had been thrust eastward over vertical 
Mesozoic sandstones, shales and coals on the Turtle Mountain 
Fault. The slide mass was reported to have moved eastwards 
downmithe Edipieommraitaset of strike Joints perpendicular to 


bedding. 


A geological map of the area was made in 1911 and 1912 
by W. W. Leach as part of the Geological Survey's program 
of mapping coal deposits. It was published with a report by 
J. OD. MacKenzie (1913) who first recognized the Turtle 


Mountain anticline while mapping south of the slide mass. 


In 1931, J. A. Allan of the Department of Geology, 
University of Alberta, was commissioned by the Alberta 
Government to investigate the conditions of Turtle Mountain. 
From mapping the geology and the fissures on top of the 
mountain he realized that if there was any danger of another 
slade; it lprobably existed at the south peak where, on the 
east face of the mountain, the bedding planes dip to the 
east. Allan (1933) presented a cross-section of the geology 
of the south peak as shown in Figure 5.3. Allan's report 
was never published but a copy exists in the Geology Library 


at the University of Alberta. 


15h Be MacKay (1932) of the Geological Survey of 


Canada mapped the coal fields in the Crowsnest area and 
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while there considered the geology of Turtle Mountain. In 
his 1932 report, he gave a cross-section of the structure of 
the south peak (Figure 5.4). MacKay recognized the easterly 
dipping beds on the east face but did not recognize the 


anticline in the mountain. 


Norris (1955) has published a map of the Blairmore area 
in which he also gives a general cross-section through 


Turtle Mountain (Figure 5.5); 


The succession of rock units within the Turtle Mountain 
area has been summarized by Norris (1955). The oldest rocks 
exposed on Turtle Mountain are Mississippian and belong to 
the Banff Formation. The uppermost 400 feet of this 
formation crop out and are medium-grained, crystalline and 
Shaley limestones containing many dark gray and pale brown 
weathering chert stringers. The basal bed of the overlying 
Livingstone Formation, 25 feet thick, consists of coarse- 
grained medium to dark gray limestone with a little chert. 
Above this are about 300 feet af massive, gray, crystalline 
limestones with interbeds of argillaceous and cherty 
limestone, the Pekisko member. The overlying Turner Valley 
member is a crinoidal crystalline gray limestone. The 


uppermost beds of the Livingstone Formation, which is about 
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1,100 feet thick here, are dark gray crystalline dolomitic 


limestones. 


“the Mount Head Formation, also of Mississippian age, 
lies conformably on the Livingstone Formation. athe. 1S 
approximately 690 feet thick. The basal members consist of 
interbedded shaly, dark gray dolomites and limestones about 
iWPe tectuethick se. There sthenutollowsmis0 toes of brown-gray 
dolomites with interbeddea argillaceous calcareous dolomite 
breccias, 159 feet of gray, crystalline limestone, 100 feet 
of limestone and dolomite, limestone and breccias and green 
Shale, and finally 150 feet of dark gray to black fine 
crystalline limestone and dolomite. The Livingstone and 


Mount Head Formations together form the Rundle Group. 


It is believed that Oovetlying Palaeozoic strata were 
not involved in the slide. Die westtata scropping out 
immediately beneath the Turtle Mountain fault belong to the 
upper portions of the Jurassic Fernie Group; here, 
sandstone, siltstones and _ shales. The gherniies Group lis 
overlain by the Kootenay Formation which is Lower Cretaceous 
in age. The basal bed of the Kootendyerommativonars a dark 
gray, massive sandstone. Above this is a _ sequence of 
Shales, sandstones, and coals containing four major coal 


seams. The Kootenay Formation is about 490 feet thick. 
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From the literature on the geology of Turtle Mountain, 
it appeared that the cross-section of the north peak, which 
had been published in many textbooks, was incorrect. If a 
proper stability analysis was to be done, it would be 
necessary to establish clearly and correctly the geology of 
this mountain. With this in mind several field trips were 
conducted in the area with two main objectives. The first 
was to determine the geological structure of the area and 
the other was to establish which fabric elements of the rock 
mass were kinematically active during the slide. Figure 5.6 
shows an airphoto of the slide area. This airphoto and the 
geological map shown in Figure 5.14 should be viewed in 
cOnjunceron Swithe the “plan, | shown ineariqurel san, fofetne 


geological traverses run in the Turtle Mountain area. 


The first traverse was run up the east face of the 
mountain south of the southern margin of the slide. Here 
are many trails cut for access to the coal workings in the 
Kootenay Formation. These trails were followed by vehicle 


as far as possible before accent on foot was begun. 


In the first exposures of the Palaeozoic limestone 
(thrust over the Mesozoic coal-bearing formations) above the 
trace of the Turtle Mountain Fault, the beds dip 74 degrees 


eastwards. Towards the crest of the mountain there is a 
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clear view southward along the Segzonal  straukesot Waiiicrest 
Mountain (Figure 5.8). In this figure, the easterly dipping 
beds which form the east limb of the Turtle Mountain 


anticline can be seen. 


Approximately 500 feet below the mountain crest on the 
eastern face beyond the southern flank of the slide, the 
narrow hinge zone of the anticline was found. tft coincided 


with the anticlinal crest. 


Views from the crest of the mountain showed that the 
west face of the mountain is formed by bedding dipping at 
high angles to the west. Looking north from the crest 
towards the southern portion of the scarp of the slide shows 
that, at the crown, the dip of the strata is less than 30 
degrees to the west. The scarp trace viewed from the south 
peak is concave to the east (Figure 5.9). Figure 5.9 shows 
the scarp to be almost vertical. It also shows the gaping 
cracks in the area between the peaks. These cracks are 
generally parallel to the strike of the beds and often gape 
three feet or more. There is also some relative vertical 
movement between large blocks outlined by this crack set and 


a set of cracks parallel to the dip of the beds. 


Allan, in his survey in 1933, measured the width of a 
number of these cracks. Some of his Survey stations were 


reoccupied but at none of them had there been perceptible 
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movements in the 39-year interval. Table 5.1 gives the 
measured widths of some of Allan's stations and those 


obtained at the same stations in the summer of 1971. 


The disturbed rock mass at the top of Turtle Mountain 
is not suitable for a detailed survey of the rock fabric, 
but the main fabric elements are evident. There are two 
joint sets developed perpendicular to bedding; one is 
parallel to the strike of the Weds and the other parallel to 


their dip. 


A second traverse was run up the northern edge of the 
slide debris to the area referred to as the "hoodoo" 
weathered area. Here a minor thrust exists above the major 
thrust (Turtle Mountain Fault). Overturned Livingstone 


Formation lies between the thrust faults. 


The” Bantfovwrormation’®above the minor’ thrust is highly 
ands coinplexl ye "detorned® @ast (shown = inevFagure!s 5.410. The 
anticlinal hinge trace here is only 330 feet or so above the 
minor thrust. tn this area, the east limb of the fold dips 


about 50 degrees to the east. 


The orientation of the minor thrust plane itself is 
difficult to determine from the limited exposure. However, 
because it terminates against the Turtle Mountain Fault as 


this flattens up-dip, the minor thrust is probably flat: 
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lying. Beneath the minor thrust at the north margin of the 
slide, the Livingstone strata are thrown into open recumbent 


folds. 


A third traverse was run up the southern margin of the 
slide to a point well above the major thrust fault. Here 
the minor thrust was not found. Above the Turtle Mountain 
Fault, Mount Head beds, mapped earlier as dipping westward 
ayes 0 degrees (Norris, 1955), are in fact overturned and can 
be traced upwards into beds dipping eastwards at high 
angles. These beds show a well-developed hoGi2on tap 
fracture cleavage that may represent the earlier stages of 
the minor thrust. Looking further up the slope, bedding 
appears to dip eastwards below the summit of the south peak 


but bedding planes do not daylight on the eastern slopes. 


An oblique airphoto! of Turtle Mountain viewed from the 
northeastersy shown. in, Figure, 5.11. The easterly dipping 
beds can be distinctly seen along the southern margin of the 
slide. In fact, when viewing the mountain from this 
advantageous angle, even the crest of the Turtle Mountain 
Anticline is readily visible in the south peak. It can also 
be seen in Hillcrest Mountain beyond Turtle Mountain. This 
photo shows the sharp crown and westerly dipping beds on the 
west face. 


1 This photograph was obtained from Chevron Standard Ltd. 
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An area which is of particular interest outside the 
slide area is the small abandoned limestone quarry at the 
northerly end of Turtle Mountain. Me schaseequarry the 
importance of bedding in controlling the form of the slope 
becomes obvious. On the west dipping limb of the anticline, 
there are some overhanging slopes up to 35 feet high which 
dip at 70 degrees. Striations can be clearly seen on the 
bedding surfaces (Figure 5.12). The striations parallel the 
dip of beds and are perpendicular to the fold axis 
Suggesting they were formed during flexural-slip folding of 
the strata. Jointing in the quarry is similar to that on 
the mountain crest; two joint sets are perpendicular to one 
another and to bedding, one set is parallel tor the®idip Once 
the beds, and the other is parallel to the strike. The same 
dominant joint sets are also evident in the larae blocks 


found in the slide debris. 


No Mesozoic rocks were found in the slide debris. From 
this it was concluded that the toe of the slide cropped out 
at or above the trace of the Turtle Mountain Fault. Many of 
the blocks within the fide debris have surfaces with 
Striations similar to those on the quarry wall. This would 
Suggest that the flexural-slip surfaces were also present 
within the rock mass which fell away. Figure 5.13 shows a 
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Figure 5.14 shows a new geological map of the slide 
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area. It is based on a map published by Norris and the 
additional detailed mapping done during the course of this 
thesis. The new map differs only in detail from the one 
presented by Norris. The additional mapping, however, 
provided information on the smaller, mesoscopic 
discontinuities within the slide area and allowed the 
identification of the kinematically active discontinuities. 
It also aided in placing the position of the minor thrust 
fault above the Turtle Mountain Fault more precisely and 
resulted in more information on the orientation of bedding 


on the east limb of the Turtle Mountain anticline. 


From the new geological map and the additional 
information obtained from the detailed mapping, it was 
possible to draw new cross-sections of the slide area. 
These are shown in Figures 5.15 to 5.17. Their location 


with respect to the slide is shown in Figure 5.14. 


In order to draw the general shape of the anticline, a 
profile was taken along a section parallel to the actual 
sections but just to the south of the slide margin. This is 
an area with a fair amount of data on the orientations of 
outcrops (Figure 5.14). Each one of the dip and strike 
readings in the vicinity of this section was projected on to 
the section along the strike. The dip readings were 
converted™? to the apparent’ dip along the Section. Theil shape 


of the fold was then drawn according to Busk's method of 
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constructing folds (Busk, 1929; Bablelaings, TS72) : 


This general Shape was then projected onto the 
sections. The positions of the folded strata were based on 
the traces of the formation boundaries, on the thickness of 


the various rock units and the position of the fold axis. 


The profiles of the mountain before the slide were 
drawn on the basis of a photo of the mountain before the 
Slide and the descriptions given by residents of the area 
familiar with the mountain before the slide (McConnell and 
Brock, 1904). It was reported that the portion which fell 
away was at a higher altitude than either the north or south 
peaks and that it protruded to the east farther than any 
part which remained. On the basis of this,*a line; slightly 
concave to the west, was drawn from the intersection of the 
contours and the slide{ (margineron cone peésidé wto the 
corresponding point on the other side. These lines were 
taken to be the contours before the slide. The pre-slide 
profiles were then drawn with these estimated contours as a 


guide. 


5.5 DISCUSSION 
The re-examination of the geology of Turtle Mountain 
has established that the Frank Slide is one in Which the 


major portion of the slip surface follows the bedding 
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planes. It has also been shown that at the toe of the 
slide, the surface of rupture follows a flat-lying minor 
thrust fault located above the Turtle Mountain Fault. The 
structure to the south of the slide mass, along the trend of 
the Blairmore Range, is similar to the structure at the 
slide except that the minor thrust is not present. = Ouest 
would seem that the presence of the minor thrust was a 
crucial factor in the stability of the slide mass. These 
observations are indeed in great contrast to the thinking 
that the slide occurred entirely across bedding which has 
prevailed in much of the geotechnical literature until the 


present. 


The flexural slip folding at Turtle Mountain is a style 
of deformation associated with thrust folding in many areas 
in the Canadian Rockies (Dahlstrom, 1970). In places where 
this folding is associated with the critical position of a 
thrust fault, there is a reasonable possibility that another 
Frenk-“Siide “night occur. Indeed, there is evidence that 
other rock slides of the same magnitude have occurred 


already (Dishaw, 1967). 
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WIDTH OF CRACKS ON TURTLE MOUNTAIN 


Allan's Values Measured 
Data Sept. 20, 1971 Elevation 
Station (Feet) (Feet) (Feet) 
2 ¥3765 1263 13410 
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Figure 5.2 Cross-section through the north peak of Turtle Mountain 
(after Daly, et al., 1912) 
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Figure 5.6 Airphoto of the Frank Slide area 
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Figure 5.7 A plan showing the general areas of the geo- 
logical traverses run 
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Figure 5.8 View of north face of Hillcrest Moun- 
tain from Turtle Mountain showing 
the Turtle Mountain Anticline in 
profile 


Figure 5.9 The crest of Turtle Mountain 
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Figure 5.10 The Banff Formation in the "Hoodoo" 
weathered area 
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Figure 5.12 Striationseon bedding surfaces 1n an 
abandoned limestone quarry at the 
north end of Turtle Mountain 


Figure 5.13 A rock within the slide debris show- 
ing a sheared surface 
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Rocky Mountain Formation 
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14 Geological map of Turtle Mounta 
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CHAPTER Va 


THE MECHANICS OF* 7H 8) FRANK) SLIDE 


6.1 INTRODUCTION 
This chapter presents the studies of two different 


aspects of the mechanics of the Frank Slide. 


The first aspect is the stability analysis of the 
slide. Now that the geological setting has been correctly 
established and the shear strength of typical 

iscontinuities has been measured, it is possible to perform 


the back-analysis. 


The second aspect is the effect of the mining on the 
INLVtiavzon of Vthe®slide? Up" "to? now, “the effect ** of” *the 
mining has been largely speculative. This question is 
considered in this chapter on a more mathematical basis by 
examining the change in stresses along the slip-surface 
which resulted from the mine opening. From the change in 
stresses, it was possible to determine the effect on the 


stability of the mountain. 
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A limit equilibrium stability analysis was performed on 
the three new cross-sections, presented in the previous 
chapter, uSing the Morgenstern-Price Method (Morgenstern and 
Prace;7**1965)?. The analysis was done for three cases. The 
first was to determine the friction angle, $, required for a 
safety factor of unity. In the second case, it was assumed 
that the shear strength along the entire slip surface could 
be represented by the peak shear strength parameters 
Obtained Ttonm the  iVexural-slip® surfaces: This was done 
because this type of discontinuity is the most prominent 
along the slip surface. The third case was the same as the 
second except that the shear strength parameters obtained on 
the joint surfaces were used along the failure surface where 
it goes across bedding. The cases were analyzed with the 
assumptions of a rock density ‘of 160 pef and zero water 
pressures. The zero water pressure assumption is based on 
field observation of no surface seepage above the major 
thrust fault. Gince the slide occurred) in April, “the only 
possibilty” of*®janyo vwater*opressure st Qwould ss thave heen from 
melting snow. It is doubtful, however, if there would have 
been enough meltwater to have a significant effect on the 
overall stability. Generallyaethen conditzonsse@on SLopngos 


Turtle Mountain are unfavorable for a large accumulation of 
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snow. Since the mountain is relatively devoid of trees, the 
persistent west Pas through the Crowsnest Pass area would 
prevent a tlarge accumulation of snow. Furthermore, since 
the strata at the crest of the mountain dip away from the 
scarp, the meltwater could have accumulated only from the 
relatively small area between the pre-slide scarp and the 
present scarp. Also if there was any meltwater it would 
only have been active near the crest of the mountain and not 
in the lower regions of the failure surface where the mass 


obtained the majority of its resistance to sliding. 


The results of the limit equilibrium analysis show that 
the average @ required for a factor of safety of 1.C is 
37.59. For cases two and three, where the analyses are 
based on the shear test results, the average F.S. is 0.86. 
The particular results for each case and cross-section are 
presented in Table 6.1. The average computed F.S. of 0.86 
differs by omly 14% from the defined F.S.0f 1.0. . Chis vic ia 
pleasing result in the light of the discussion given in 
Chapter 2 where it was shown that the difference can often 


be several hundred percent. 


Another way to inspect the agreement between the 
analysis and the defined F.S. of 1.0 is to calculate the 
average normal and shearing stresses on the failure surface 
andiacthen s plot them in conjunction with the shear strength 


envelope. If these average stresses fall above the failure 
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envelope, Stherer.s? wouldsoberlessPrhanione.i8 Tha spwase done 


for the Folkestone Warren Slide by Hutchinson (1969). 


The foregoing procedure was used on the three new 
cross~sections “of “thet Frank "Slide? 9The average” stresses 
were calculated at the base of each slice on the basis of 
knowing the weight and side forces of each slice. The side 
forces were taken as calculated by the Morgenstern-Price 
program ata F.S. of 1.0. The stresses at the base of each 
slice were then averaged over the whole slip surface. The 
results are shown in connection with the peak strength 
envelope in Figure 6.1. The average stresses fall just 
above the strength envelope indicating that the F.S. is less 
thane 1..0), This is in agreement with the computed safety 


factors shown in Table 6.1. 


Presenting the stability analysis this way allows the 
results to be viewed in conjunction with the scatter in 
shear test data. As can be seen in Figure 6.1, the average 
stresses for the Frank Slide failure surfaces fall just 
above the strength envelopes. They are close enough to be 
within the 95% confidence level of the mean of the peak 
flexural-slip shear test data points. Furthermore, the 
average stresses generally fall within the scatter of the 
test data. This is considered an excellent result for the 


analysis of rock slopes. 
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The assumption of zero water pressures is supported by 
the relatively good agreement between the F.S. obtained 
using the shear test results and a F.S. of 1.0. The fact 
that this variable is not included in the analysis is 
probably one of the reasons for the good agreement. Many 
times in a back-analysis, assumptions have to be made about 
the water pressures at failure. If these assumptions are 
incorrect, it is not possible to tell if the shear tests are 
Depresentativciotethe field iconditions.ef Thistieimnot tov gay 
that any small increase in the water pressure due to 
precipitation could not have acted as a triggering mechanisn 
because, as has already been demonstrated, the stability of 


the entire mass was in a very critical state. 


Generally, the F.S. in rock slope analysis is very 
sensitive to a change in the cohesion, C. This is not so in 
the analysis of the Frank Slide as shown in Figure 6.2, the 
reason for this being the height of the slope. For a given 
slope geometry anda given # and water pressure condition, 
the F.S. will depend only on the C/yH ratio (Bishop and 
Morgenstern, 1960). H is the height of the slope and y is 
the bulk density. In the case of a small H, the change in 
C/y H will be much larger for a given change in C than at a 
high» +H. Since the F.S. is a function of C/Y H, the same 
degree of change applies to the F.S. for a given change in 
ee At the Frank Slide where the height of the slope is 


approximately 2,000 feet, the F.S. is therefore insensitive 
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to a change in C. This is probably another reason for the 
relatively good agreement between the computed F.S. and 


Wry ¢ 


The stability analysis showed that the average ) 
Tequired Olea tr. co. Of 1eO,) "Withee 2zeLe Coresion, Bers 3750 
degrees. This is very close to the basic friction angle of 
37.29 measured on the surfaces lapped with the #45780 grit. 
This type of agreement has been observed on other slides as 
discussed in Chapter 2. On the basis of these observations, 
it is the writer's opinion that a fairly good first estimate 
of the degree of stability can be made by using the t, as 
measured on a flat surface which is sandblasted or lapped 
with a fairly coarse grit and with C equal to zero. siisiibis 
would provide a lower bound to the problem and be on the 


safe side. 


It is interesting to note that ¢ required for F.S. of 
Alpalls ali Sho ss ga vope 6, + i but simply oe Furthermore, the shear 
strength as determined on the flexural-slip surfaces has 
been used in the stability analysis with good success 
without any correction for the field roughness. This would 
seem to indicate that the effective i at the field scale is 
essentially zero. Again, to assume the i to be zero would 


present a lower bound solution on the safe Side. 
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65 3) THE) BFFECT OF, MININ 

MOUNTAIN 

At the time of the Frank Slide, the Canadian American 
Coal Company was mining a coal seam at the base of Turtle 
Mountain which was approximately 15 feet wide and dipped 
approximately 85 degrees to the west. This coal seam is 
identifiable in Figure 5.6 where the entrance to the mine is 
at the junction of the north margin of the slide and the 
river .f.On thee south edge) of..jthe:.slide,,: subsidence _is 
Clearly visible marking the seam at this point. Mining had 
reached the south edge of the slide debris when the slide 


occurred. 


The coal was mined in a series of chambers about 130 
feet in length and and from 250 to 400 feet high. These 
chambers were separated by pillars 30 feet wide which 
contained manways and timber chutes. These large chambers 
had reached the south edge of the slide when the disaster 
occurred. In fact the edges on both sides of the slide 
corresponded very closely With otthejpadaimite,rofe thes, bag 


chambers in the mine. 


For about six months prior to the slide there had _ been 
a general squeezing in of the tunnel walls (McConnell and 
Brock, 1904). The coal would mine itself by breaking off 


the hanging wall. The miners noticed slight movements at 
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times which were associated with shocks which felt like 
tremors. GangwayS and manways often had to be retimbered 
and were even condemned due to the impossibility of keeping 


them timbered (Daly, et al., 1912). 


Exactly how much movement occurred in the mine is 
difficult to determine. Daly, et al. (1912) have, from the 
amount of coal removed, estimated that if the entire void 
left by the coal was squeezed in, it would represent about 
four feet of movement. Although this could happen with 
time, it was not the case just after the slide occurred 
because these chambers were mined again after the slide, and 
witnesses also reported that the mine escaped with little 
damage. The’ testimony’ sof the miners at™the®  timevoft “the 
disaster contains nothing that would indicate that there was 
a sudden movement within the mine walls just before the 


slide. 


TENLIS, SORAECOULTSEFRinportant to Ssstablish what effect 
theurmining had cone dthes initiation” of the slide Lf 
recommendations are to be given on similar future projects. 
This was investigated by determining the stress changes 
resulting along the failure surface using the Finite Flement 


Method. 


The stress analysis was done with a Finite Flement 


computer program based on a simple constant strain 
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triangular element. The mountain was divided into a finite 
element gridgastshownein Figure 6.3 (Section ¢-C"'). The 
initial insitu stresses were determined using the "gravity 
turn on" method (Goodman and Duncan, 1968b) where K,, the 
coefficient of earth pressure, was set equal to V/1-v (v= 
POlsson'S La v10) @ SOniy the Condition for Kj) = Us is" was 
analyzed. From the results of this analysis, the initial 
insitu normal and shearing stresses were calculated alorg 


the failure surface (Figure 6.4). 


The mining at the base of the mountain resulted in the 
walls of the mine becoming a stress-free boundary. The 
effect of this in the stress distribution can be simulated 
by applying stresses to the mine boundary equal and opposite 
to the initial insitu stresses. With these reverse stresses 
applied to the mine boundary and the mass taken as 
weightless, the results from the analysis gives the change 
in stresses throughout the mass due to the mine opening. 
The finite element mesh with the mine open is Shown in 


Figure 6.5. The change in stresses is shown in Figure 6.6. 


The effect of these stress changes can be evaluated by 
defining 4 pseudo safety factor, SF , as the resisting 
forces along the failure surface divided by the shearing 
forces. If we assume C = 0, the resisting forces would be 
tan $5, as integrated over the whole failure surface. Ono 


is the normal stress over an incremental distance, ds. The 
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Shear force integrated over the failure surface would be 


T, ds. In equation form 


b 
f tan PPG ids 
a O 


SFo = 


@eeeeceeesevseeseeeee ee 6. | 


Mee 
a 
lk 


where a to b is the distance along the failure surface. 


After mining, the stresses along the failure surface 
re 
would be On, = ( 7 Oeeeand) Cra +47, where a6, is the 


change in normal stress and47Z is the change in shearing 


stress. The safety factor after mining would be 


b 
tan fg iP iGe nara) cls 
a 


of = 6.2 
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b 
f (7, +47) ds 
a 


The change in safety factor then is 


From the distribution curves (Figures 6.4 and 6.6), we 
can integrate using Simpson's rule. The results of these 


integrations are as follows: 


ie ds = -100.2 x 103 Kips/ft 


a opal © wree 


Pid cv awen wees sae eetuns nigel eee M of B py Powers 
wah *. Veen bag 


5 
esyiiecs 
eoetive ozplish eft pnolts somegenb odt ei o eo? 8 stedv 
» €s3tee 
Ve 


soptis e1list eds profs eseze7se edt ,pninia ett 


eit at ,0s eTedw a+ 5) = a bas aeat on = .av od Bivow 
patragie at opasd> off at Yo has ses1#e Lemzon ab epasdo 
aids 
od bivow patcte reste wos Det at oer eee tte 
j rae ey 
fo,3 : i 
2b (,04+ > \ ost 
a 
< 8 pa aq Y 
C.o sseveeeeaeveere = 2 
: abe. 
ab ($44.3) 
8 , "41¢e4N 


ef aed? toss? yistsa mi epasd> oft 


je = ge s Of 
: - 7=44 
ew (8-3 Bas 8.3 eeuuees) zevIED nottudizsetb od? a 


ezeant to etineea oiT .oIDz a"noaqni2 priev ssaviesal 18> 
: AL aLSee 
swt 7h hs ve enoltts1petal 

ake >. i 


Preait €or x S.00T- = 8b, s\ | 
Ss ee cay ond: othr 


é¢ ops * \qatee cle eb 


‘a 


7 ae aa 


156 


Si 
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69.9 x 103 Kips/ft 


1,39 °x 103 KipsyeEt 


m@ Ss 
> 
a 
Qu 
(67) 
tt 


i fen 00 10s Rane ee 
a 
In numerical form then 


SPS =-can 0 100.2 
69.9 


an! 9° (1.433) 


SF = tan @ (-100.2 + 1.36) 
( 69.9 = 0.20) 


= tan @ (1.418) 
The percent change in SF then is 
tan # (1.433 - 1.418) x 100 ¥ 1% 
tan (1.43) 
With this method of analysis, it is also possible to obtain 
@ for "a factor of Safety of unity. If we set oF equa. 1o 
1.0, we can solve for @ from 
SF, = tan $ (1.433) 
WLDNSESES eqtal 8t041)0 


O8=Gtan st 


(Bees) 
(1. 433) 
Y 350 


This is in good agreement with the values obtained from the 


Lins Sequijabrium stabvlity analyses (Table 6.7). 
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fixed boundaries, for a two-dimensional anaylsis in a 
homogeneous isotropic body is not affected by the magnitude 
of the elastic modulus, E. Thus for the above calculations 
of the stress changes, the results are not affected by the 
choice of E. To assume the entire mountain to behave as a 
homogeneous isotropic body is a simplification of the 
behaviour of the material but does allow an estimate to be 


made of the effect of the mining. 


The deformations resulting from the mine opening, 
however, are directly dependent on the magnitude of E. The 
inward deformation at the top©® of @tthe! imine on © othe 
mountainside (Figure 6.7) is approximately three inches if 
we assume E = 1,000,009 psi. Decreasing the F value to 
100,000 psi increases the deformation at the same point ten- 


fold or to approximately 2.5 feet. 


The general over-all style of deformation of the entire 
mountain is shown in Figure 6.8. Here we see that the top 
of the mountain in the region of the failure surface moved 
more or less as a unit. This type of movement accounts for 
the low change in shear and normal stresses along the 


failure surface. 


Generally, this stress analysis does not indicate that 
the mining was a factor in significantly decreasing the 


stability of Turtle Mountain. FLOM selina equilibrium 
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analysis, however, it has been shown that the stability of 
TurtlesMounrain, Wasi in) a. very . critical ostate,,. heforas the 
slide sand with. this an mind, the 1%, decrease anustabi lity 


due to the mining may have been enough to trigger the slide. 


6.4 DISCUSSION 

The stability analysis together with the shear strength 
results have shown the Pees unstable state of MTurtle 
Mountain before the slide occurred. What finally triggered 
Thess Wide. Lsndat ti cult bo. say. Factors which have been 
cited as possible causes are the ice wedging in the cracks 
on top of the mountain, earthquake tremors, increased pore 
pressures and mining at the base of the mountain (McConnell 
and Brock, 1904). All these factors, although small, cannot 
be ruled out completely as being responsible for snapping 
the last bond by which the mass was held up. The primary 
cause, however, was undoubtedly in the form and structure of 
Turtle Mountain. I+ might be argued that this is not new 
Since it is the same conclusion that was expressed by 
McConnell and Brock in 1904. What is new, though, is that 
now it has been substantiated by a stability analysis, and 
the new interpretation of the geology brings more logic to 
why the structure waS so important. The added advantage of 
having shown this from a stability analysis is that the 
analysis can now be used with a greater degree of confidence 


to estimate the stability in similar situations. 
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One of the main items of concern throughout this thesis 
has been the reliability of the stability analysis of rock 
slopes based on shear test data. The analysis of the Frank 
Slide indicates that for a slope where the slip-surface 
follows some through-going discontinuity, where the material 
is a hard unweathered rock, and where the C/YH ratio is 
small, a fairly reliable estimate can be given of the 
F.S. on the basis of laboratory shear tests. This 
statement, however, assumes that the geological and 
geometric setting are known in detail and, furthermore, that 
the water pressures at the time of failure have been 


correctly established. 


In recent years, there has been an increased interest 
in developing recreational facilities along the eastern 
Slopes of the Canadian Rockies. Along these slopes, there 
is a possibility of encountering a geological setting 
Similar to that of Turtle Mountain. The understanding 
gained about this type of situation from the Turtle Mountain 
investigation would be very valuable in making 


recommendations on areas to avoid for planned development. 
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TABLE 6.1 


RESULTS OF STABILITY ANALYSES 


ON FRANK SLIDE 


i 


Case l Case 2 Case 3 


F.S. for #=28° 


@ required for Cc=4600 psf and 
F.S. equal to F.S. for f=28° g=329, C=3600 psf 
Section thers) C=4600 psf across bedding 
A-A' 39.7 0.832 0.837 
B-B! B52 00 0.890 0.896 
CeCe 37.4 0.853 0.863 


NOTE: (1) The average for Case 1 is 37.5 


(2) The average for Cases 2 and 3 is 0.86 
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Figure 6.1 The average normal and shear stresses in relation 
to the failure envelopes 
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SHEAR STRESS 


KIPS 7 FT ° 


NORMAL STRESS 


MST 80 I6O 240 


DISTANCE ALONG’ SLIP SURFACE ( FT. X 10) 


Figure 6.4 The initial stress distributions along the 
slip surface 
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Figure 6.6 The changes in stress distributions 
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Figure 6.7 The inward movement of the mine walls 
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CONCLUDING REMARKS 


In this thesis, an examination was made of the 
reliability of rock slope stability analysis by inspecting 
humerous case histories and by performing a detailed back- 
analysis of the Frank Slide. Thais “study? ocfieetai Vediierock 
slopes showed that in the majority of cases, the shear 
strength parameters required for a safety factor of 1.0 are 
# equal to 30 + 5 degrees and C equal to zero. Furthermore, 
the case histories revealed that the F.S. based on shear 
tests can differ from the defined safety factor of 1.0 
Within a range of only a small percentage to several hundred 
percent. For this reason, rock slopes are often designed on 
the basis of past experience in conjunction with a 
consideration of local geology and not so much on stability 


analysis. 


ine *connectton swith » they back-anallysis Jot * «hel Prank 
Slide, shear tests were performed on natural discontinuities 
and artificially prepared limestone surfaces. The results 
of these tests yield not only the necessary strength 
parameters for the stability analysis but also shed some new 
light on the shear strength behaviour of rock 
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dL Scontinuuwteves. 


An investigation was conducted into the possibility of 
mathematically characterizing the microscopic surface 
roughness of prepared rock surfaces and then correlating 
these with the friction angle. The results of this work 
showed that it is possible to establish this type of a 
correlation; however, the work involved in calculating the 
roughness parameter makes this an impractical method for 
determining 9@. TLS) application. to the, Viiteldiwscale as 


hindered by the inability of measuring the field roughness. 


The shear tests on the natural limestone 
discontinuities revealed that the shear strength is not 
affected by the size of the sample if the scale of roughness 


does not change with the size of the sample. 


The ultimate shear strength (shear strength after 
larger shearing displacements) was measured on the natural 
discontinuities by performing a reversal shear test whereby 
the geometric component of the shear strength was 
eliminated. The results showed that the ultimate strength 
is greatly affected by the initial roughness of the 
discontinuity and the type of surface damage which occurs 
during the shear tests. Since the ultimate strength is a 
function of these conditions, it differs characteristically 


from the residual strength of clays. 
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The re-examination of the geology of Turtle Mountain 
proved to be most enlightening. A cross-section showing the 
bedding dipping into the slope and the failure surface going 
across bedding has appeared in many textbooks and 
publications.® "This cross-section is incorrect The) strata 
within Turtle Mountain are in reality folded so that on the 
southern margin of the slide, the bedding dips in the 
easterly direction out of the slope. With the strata folded 
in this manner, the failure surface only went across bedding 
at the crest of the slope. Further down it follows the 
flexural-slip surfaces parallel to bedding and then along a 


Minor thrust fault at the toe. 


The average “safety “factor off the 7FPrank “Slide as 
computed on the basis of the shear tests differs, on the 
unsafe side, by only 14% from the defined safety factor of 
120. *For rock slope “anallysis, “this "is"considered ‘to be very 
good. “This indicates “that = for®™a “potentials slidelgin Ga 
similar geological and geometrical setting, it is possible 
£6" "predict the’ degree stor ‘stability with considerable 
reliability from shear tests performed ony indt ural: 


discontinuities. 


The study of the case histories has shown that with 
C™="Qe"8the ¢ required for a F.S.s of 1.0 is very similar to 


the basic friction angle, #,. This suggests that a good 
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The stress analysis of Turtle Mountain demonstrated 
that the mining at the base of the mountain had a small 
effect ane decreasing Sthe Sstability BRorlythesal rocthomass. 
However, Since the potential sliding mass was already in a 
critically unstable position, it is not unreasonable that 


the mining triggered the slide. 


With regard to future research in the area of rock 
Slope stability, there is still a need for further detailed 
case history studies. The reason for this is two-fold. 
Firstly, a back-analysis is like analyzing a large shear 
test and, therefore, is a means of investigating the shear 
strength behaviour of rock discontinuities at the field 
scale. For example, the study of failed slopes in this 


thesis seems to indicate that if there is a problem with 


instability, the effective i at the field scale is 
essentially zero. Only the analysis of failed slopes can 
provide this type of information. The second reason for 


further case studies is that they would clarify the 
reliability of analysis based on shear tests; i.e., they 
would either confirm the present position or provide more 
definite bounds on the reliability for various geological 


settings. 
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Tn the overall Systeme of “stability )analysisc,, othe 
problem is not so much the type of analysis as it is the 
establishment of the field shear strength parameters. 
Therefore, it is suggested that future research efforts be 
directed at understanding the field shear strength behaviour 


rather than the development of any new analytical methods. 


In the design of open pit slopes, a change in even one 
degree in the overall slope angle can represent either a 
cost or saving of millions of dollars. Ideally, it would be 
desirable to design the slopes within this low tolerance and 
be certain that there would be no major failures. This, 
however, SS “an unrealistic Vobjedtive, Jandi) ip )ais not 
suggested that further studies of failed slopes will provide 
the necessary confidence in stability analyses to design at 
this low tolerance. The nature of the material within rock 
slopess isi. too, complex Sim behaviour ) for? thisiwiita, bea 


realistic objective. 
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APPENDIX A 


OPTICAL POUGHNESS PLOTTER 


A.1 INTRODUCTION 


At the University of Alberta hydraulics research 
laboratory, a device has been built for plotting the flume 
bed profiles under water and yet not disturb the river bed 
sand. A device based on the same principles as the river 
bed plotter was constructed for the purpose of plotting 
roughness profiles: oOferock: discontinuities. )This¥Appendix 
Gives! ae description ® of “the: *ceonstruction;, soperation “and 


performance of the roughness plotter. 


Basically, the device operates on a light source which 
remains a constant distance from the surface. Inside the 
light probe (Figure A.1), a small light bulb is located 
whose light is transmitted along optical fibres and 
projected Onto BtHeEMMrOCK. Suniace. Mixedisinigwathy whe 
transmitting optical fibres are an equal number of fibres 
which carry the reflection from the rock surface to a photo 
resistance cell. The voltage drop across the photo cell 
depends on the intensity of the light directed at the cell. 
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The intensity of the light reflection to the fibres will 
depend on how close the tip of the fibres are to the rock 
surface. In this way the voltage drop across the photo cell 
is a function of the distance between the light probe and 


the rock surface. 


Initially, a current is passed through the photo cell, 
and this is taken as a reference point. Then as the light 
probe traverses a surface and the tip moves closer or 
farther away from the surface, a change in potential occurs 
across the photo cell. This change in potential from the 
initial reference point activates a small D.C. electric 
motor which moves the probe height until it again reaches 
the initial distance from the surface. In this way, the 
electric motor lifts or lowers the probe in such a way that 
it is always a constant distance from the rock surface. 
Basically then we have a system in equilibrium designed so 
that if any deviations occur from the point of equilibriun, 
the system will automatically adjust itself to return to the 


equilibrium position. 


A.3 CONSTRUCTION 


The entire assembly of the light probe is shown by the 
schematic in Figure A.1. Thee Helectric, BNOtOL TSseeawt 12, 
DC MOtLOL. The purpose of the bellows is to rectify any 
misalignment which might occur between the large pulley and 


the MOtor. The potentiometer measures the vertical 


Iles otody od? veotos 407 opatioy 
bunt sdorq daptt sit nisaihciost AB! 


iiss otety sult dpwolrd? boeesy at spill 6 seiteliar on: 

Fdpll odd es nodt .«tatog ssacugieris. 2s noast ak eld? ba ae a ‘ 
"0 tenolf>o aevose git ear 5ns eps Ive ss eoatevemt edomq ’ | 18 
aqwnoo Isitnetoq at epssado a ,e0s?70e edt acm yews rodsa8> | ” 


oat wot? Isttsasoq at onder eid? shies otodq en2 a raeth aii 


sittoele .>.¢ dfsae & eetsviios Pnkog eonetete7 isrenae ; 
aedosor nisps Fi Ltiau tdpied sdetq sf#2 eves doidw aotom rn | 
edt «ysw etAt ot sont rps id wom aoa steke isisigt edt . ; 
tadt qe 6 ive at odorg ont erewo! wo af2tL x9t0" a | 
,soetz¢2 Avot of? aot? sonsteth sasten0o & eqawie at 2 i 


o2 Banpieeb guiadttivups at neteq2 5 ovsd sv asit yifsolesd A: oa 
<muitdifivps to tnioq edt ao72 209900 anoitsived yas at edt A, a 
sat oF amwsst oF tieati savths yilsottsaotus Iliv graded edt ok 

.1ottteog avizdtilups eh 


. >< Zz 


vSt s 2 azotom st3topie edt. ree ezpit mi. otsenetoe A 
yne Y2itosa oF af ewoLied od? 30 s80gTHG ed? .totc 


185 


displacement. 


The light probe is mounted on a carriage which in turn 
TSiepart woteythe: holding frame (Figure Aw2). This carriage 
runs along the one-inch square rails under the power of an 
electric motor. The potentiometer under the motor measures 
the horizontal displacement. The purpose of the twelve 
roller bearings is to allow the carriage to run along the 
rails regardless of the position in which the holding frame 
is mounted. In order to accommodate the transportation of 
the apparatus, the frame was constructed so that it could be 
taken apart and reassembled with relative ease by using 


Slotted bolt holes and wing nuts. 


The electric circuit, the power supply and the control 
panel are all contained in a special box. The power supply 
consists of six 6-volt dry cell batteries except for the 


light source which uses four flashlight (1.5v) batteries. 


The entire apparatus is shown in a laboratory set-up in 
Figure A.3, and a photo of a field set-up is given in Figure 
AS. -S The field’ set-up! "shows ‘howlit cane be’ mounted on “an 


overhanging rock discontinuity. 


Recording the output from the potentiometers proved to 
be the most difficult as far as field work was concerned; 
mainly because no suitable analog recording device could be 
obtained that did not require a 110v power supply. 


Consideration was given to either an xX-Y plotter or a 
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magnetic tape recorder. The magnetic recorder was chosen 
over the plotter because it is more rugged. The recorder 
was a Hewlett-Packard Model 3960, 4-track, magnetic tape 
recorder. A small 100-watt generator with a gasoline engine 


was used as the power supply in the field. 


After returning to the lab the signal from the tape was 
fed into the X-Y plotter yielding roughness profiles in 


graphical form. 


The icalubration factonss for hewpotentiometers are as 
follows: 
vertical; 3.0125 an/V 


horizontal; 8 85a AN, 


A.4 OPERATING PROCEDURE 


Figure A.4 gives a view of the control panel for the 
operation of the optical roughness plotter. The panel is 
divided by a dotted line and only the left side of the 
control panel is concerned with the operation of the optical 
plotter. Generally, the following steps should be observed: 

(1) the button on the dividing line has a white spot on 
it and should be pointing to the left of the line. 

(2) flip the two toggle switches side by side to the “on™ 
position. 

(3) the damping gauge should be set to approximately 80. 

(4) lower the light probe manually to about half an inch 


from the surface by moving the "up-down" switch to 
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the "down" position. 

the gain should be set at approximately 10 to 20. 

the probe height should be set at approximately 80. 
flip the toggle switch "manual-auto" to the automatic 
position. With this action, the probe should move 
down to the equilibrium position. If the probe moves 
away from the surface, the probe height is set too 
high. If nothing happens, increase the gain very 
Slowly. When the probe begins to chatter and 
vibrate, the gain has been set too high. 

once a stable position has been reached, the probe 
height can be adjusted to the desired height. 

when the light probe has been set, the next step is 
to traverse the rock surface. This is done by 
invoking the toggle switch, in the lower right-hand 
corner, to the “forward” or "reverse" position 
depending on where the carriage is located. There is 
also an adjustment with which the carriage proceeds. 
The carriage will continue to run along until it is 
stopped manually or until it hits the limit switch. 
To reverse the direction, place the direction switch 
in the opposite position and depress the "override" 


button until the limit switch is free. 


SULTS OF TRIAL TESTS 


Figure. A.i6° shows thes results? offvar) traverse onan 


ai _ ye aie, 
“66S Or or f svoonancng 
(08 yLevemtkorqas ts pa 
dissgosus ont OF YOduE- ae 
svat bluoda sdozq edt ,cobdss . 
eover edoaq sif 3 .aottieog it nd £4 
oot tea ef sdpted oiorg edt \s0A3 - 
yrev sisp edt seseipdi § ,aneqga cAlthe ar : 
bas vested> of antped sdeotq sit nodW eivete . 
.bptd oot tee ased sad aisp of? ,eteadiv’ 7 en. 
edotq sft .beroset asod asd woitieog sfdsde 6 eono 8) - 
.tdpted berieeb add of Setewhbs od ono Sdpfod it : 
at qofa txen edt toe osed asd sdorq tdpkt ote asi ron ar | 


yd sob az etd? .soataee door edt semevezt OF = | 


hasi-tdptt revel eft at ,Aptiwe efppot eit prizoval — 4 a 
gottisog “sarg¥eT" 2o "haswroI" oft oF .i9mt00 : 7 . : 
ak sted? .betsool ef spsiazso sdt otedv a0 pnibaeged: nase 7 
s2besootg spaitzes edt Asidw ete tnonsevt 5s as cela x a | - ; 
2i ti Litany patois net of ormitaoo Ew epsizas® oat a 7 a 
.dotive tiwkl od? attd dt Livan zo yiispaam beqqose ~ male tee 7 


dotiwa goitoeath sit sosiq ,f02In91f5 ede eateves fOr 7? 
} 
"ehtaievo" d+ easigeh bas st Bag sine bb ott ak yh 


188 


artificial surface cut from a piece of lumber. Figure A.7 
shows the results of the runs made at the set-up shown by 
Figure A.5. These profiles show the scale of roughness 


which this plotter is suitable for and able to measure. 


A.6 PERFORMANCE 


The main advantage of the optical plotter is that the 
profile record is continuous as opposed to discrete readings 
at numerous points. Also the plotter can accommodate fairly 


large vertical displacements (10 in). 


The accuracy of the measurements are controlled to some 
extent by the distance the probe tip is away from the 
surface. Increasing the distance results in the probe heing 
unable to detect the reflections properly. Also the greater 
distance results in the light being spread over a larger 
area which is undesirable as the measurement is shown as an 
average over the area covered by the light. The accuracy 
depends as well on the consistency of the reflective 
properties of the rock surface. A change from a dark to a 
lighter surface will slightly change the distance between 
the probe tip andthe surface. For this reason, the rock 


shown in Figure A.5 was painted white. 


The most serious problem in using the optical plotter 
in the field is the interference of the sunlight. The 
reflections off the surface due to the sun simply drive the 


probe away. For this reason, the measurements made in the 
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field were done late in the evening and even after dark. 


Defects in the rock surface like cracks or sharp rises 
create real havoc with this system. As the light probe 
travels over a crack, it is simply driven into the crack 
Since there is no reflection. When the probe comes upon a 
sharp) riseman height, it iwiily simply xsun up againsuy the 


abutment and the rod on which it is suspended will bend. 


Tt would appear that the only way to overcome these 
obstacles is to stop the horizontal movement and manually 
ride over them. These manoeuvers show up on the record and 
as such have to be noted and later corrected to obtain a 


true profile. 
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Figure A.1 Schematic diagram of light probe for the field 
roughness plotter 
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Figure A.2 Schematic diagram of holding frame for the 
field roughness plotter 


Figure A.3 The optical roughness plotter 
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Figure A.4 Control panel for the optical rough-— 
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Figure A.5 Roughness plotter set-up in the field 
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APPENDIX B 


SHEAR TEST DATA 


This Appendix presents the results of the shear tests 
performed in connection with the analysis of the Frank Slide 
as srelerred sto, in “Chapter "a. Table 8.1 gives the 
designation of each sample type, the area of each sample and 
the normal load at which each sample was tested. For each 
test performed the shear load versus deformation plot and 
the associated vertical deformations are presented in 


Figures Bel to 8.31. 
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TABLE B.1 SHEAR TEST DATA 


NORMAL 
AREA LOAD TYPE OF 
SERIES NUMBER (in sq) (lbs) SURFACE 
1 29.9 2,494 
2 2955 4,491 
3 30.9 7,495 
Bedding 
B 5 Sél2 556 planes 
6 3.03 1,146 
¥/ 3620 1,720 
8 BS e42 854 
1 3035 2,695 
2 31.4 4,693 
3 30.6 7,694 Flexural- 
slip sur- 
4 S0e2 10,695 faces 
5 3.03 556 Note 
5A 3.2.03 556 (5 & 5A = 
S the same 
6 Bee 854 sample but 
run in dif-— 
6A ea 854 ferent 
directions) 
t) 3.02 1,146 
TA soe 1,146 
8 Be O02 1,719 
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TABLE B.1l (continued) 


NORMAL 
AREA LOAD TYPE OF 
SERIES NUMBER (in sq) (lbs) SURFACE 
1 30.3 2,894 
2 3151 4,892 
J Joints 
3 cyeuane 1,393 
4 32! 10,880 
1 3261 3,294 
2 Sea, 5,294 
D Diamond—-cut 
3 31.0 aypeac) |) surfaces 
4 S22 in 2 ae 
1 Eis S| 3,095 
2 334 5,098 
L Lapped with 
3 S25 8,095 #45/80 grit 
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Figure B.1 Shear load and vertical deformation versus 
the horizontal deformation for sample B-l 
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Figure B.2 Shear load and vertical deformation versus 
the horizontal deformation for sample B-2 
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Figure B.3 Shear load and vertical deformation versus 
the horizontal deformation for sample B-3 
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Figure B.4 Shear load and vertical deformation versus 
the horizontal deformation for sample B-5 


0.20 


(X10-} ) 
0.10 


VERT DEF (IN) 
200.00 240.01 -0.20 -0.10 -0.00 


} 


160.00 


(LBS (X10! 


120.00 


SHEAR LOAD 
80.00 


40.00 


0.00 


Figure B.5 
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Shear load and vertical deformation versus 
the horizontal deformation for sample B-6 
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Figure B.6 Shear load and vertical deformation versus 
the horizontal deformation for sample B-7 
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Figure B.7 Shear load and vertical deformation versus 
the horizontal deformation for sample B-8 
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Figure B.8 Shear load and vertical deformation versus 
the horizontal deformation for sample S-l 
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Figure B.9 Shear load and vertical deformation versus 
the horizontal deformation for sample S-2 
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Figure B.10 Shear load and vertical deformation versus 
the horizontal deformation for sample S-3 
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Figure B.1l1 Shear load and vertical deformation versus 
the horizontal deformation for sample S-4 
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Figure B.12 
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Shear load and vertical deformation versus 
the horizontal deformation for sample S-5 
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Figure B. 13 
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Shear load and vertical deformation versus 
the horizontal deformation for sample S-5A 
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Figure B.14 Shear load and vertical deformation versus 
the horizontal deformation for sample S-6 
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Figure B.15 


213 


0.20 0.50 0.60 0.70 0.80 


0.30 0.40 
HORIZ OISP (IN) 


Shear load and vertical deformation versus 
the horizontal deformation for sample S-6A 
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Figure B.16 Shear load and vertical deformation versus 
the horizontal deformation for sample S-7 
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Figure’ B.17 
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Shear load and vertical deformation versus 
the horizontal deformation for sample S-7A 
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Figure B.18 Shear load and vertical deformation versus 
the horizontal deformation for sample S-8 
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Figure B.19 
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Shear load and vertical deformation versus 
the horizontal deformation for sample S-8A 
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Figure B.20 Shear load and vertical deformation versus 
the horizontal deformation for sample J-l 
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Figure B.21 Shear load and vertical deformation versus 
the horizontal deformation for sample J-2 
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Figure B.22 Shear load and vertical deformation versus 
the horizontal deformation for sample J-3 
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Figure B.23 Shear load and vertical deformation versus 
the horizontal deformation for sample J-4 
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Figure B.24 
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Shear load and vertical deformation versus 
the horizontal deformation for sample D-l 
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Figure B.25 Shear load and vertical deformation versus 


the horizontal deformation for sample D-2 
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Figure B.26 Shear load and vertical deformation versus 
the horizontal deformation for sample D-3 
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Shear load and vertical deformation versus 
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Figure B.28 Shea 
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Figure B.29 Shear load and vertical deformation versus 
the horizontal deformation for sample L-2 
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Figure B.3l1 Shear load and vertical deformation versus 
the horizontal deformation for sample L-4 
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SPRGl RA DENSITY COMPUTER PROGRAM 


Special techniques are required for calculating the 
Fast Fourier Transforms (FFT) of a data sequence. To 
perform these complicated calculations, a computer program 
was obtained from the Physics Department, University of 
Ribera. ee hictind = Ofmchivcr = prOgTaN el Ss  eGLVSN sin nUnis 


Appendix. 


Once the FFT are obtained, it is relatively a simple 
matter to calculate the spectral dersity estimates. The 
procedures for determining these estimates have been given 
by Bendat and Piersol (1971) and are given below in the 


context of the specific problem at hand. 


The spectral density estimates, G,, are defined as 


G(s) = 2 e h/N a Fa pte Be ee ia re ale Co) 


where h = the distance or time between amplitude readings 
N = total number of readings 
aie = k/Nh k = OUMig2 e e e v N-1 


seas the Fast Fourier Transform of the data. 


The procedural steps followed an calculating the 
spectral density estimates are as given below: 
(1) The data was truncated to a total of 2048 readings. 
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The FFT method requires that WN = o'aiwmheny P = 11, 8 
= 2048. 

The data sequence was tapered at the ends Ry ~ Using 
the cosine taper window shown in Figure C.1 
(SUBROUTINE COSTAP in the computer program). 

The ale values were calculated by first calling 
SUBROUTINE ONE R FT (2048, x), which calculates the 
real and imaginary parts of X, and stores them in the 
X atray. Then, from the New xX atray, the pe Gua le 
values were computed from the following equations: 
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A total of N/2 + 1 estimates can be obtained this 
way. 
The spectral density estimates were then calculated 
according €o equation C=1 with f = 0.0001 inches. 
The frequency corresponding with each G, estimate was 
calculated from 

f = k/N/h 
The G, estimates were then scaled by a factor of 
VW7Oe! Sedue tothe wini vial cosine tapering. 
The raw G, estimates were smoothed by using 4 
Trapezoidal Spectral Window. The smooth estimates 
are given hy 
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In the computer program, the frequency smoothing is 


done in SUBROUTINE FREQ. The Lineal 


smoothed 


frequency curve was then plotted on the Calcomp 


Plotter by calling SUBROUTINE CGPL. 


The input information required for the program 
the x array of amplitude readings taken from 
file. 

N8, the number of sets of data. 

SAM, the sample number. 

LB, the 1 value. 

NP, the number of data points to be printed. 


VB, the vertical scale on the plots. 
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all Ae O T/2 


Figure C.1 Cosine taper data window (after Bendat and 
Piersol, 1971) 
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DIMENSION X(2048),G(1045) 
LOG2N=11 

N=2**LOG2N 

READ (5,75) N8,SAM,LB,NR, VB 


75 FORMAT (15,1X,A4,215,F10.0) 


DO 100 N7=1,N8 


READ DATA FROM FILE 


DO 10 I=1,8 


10 READ (1,1) X(I) 


1 


FORMAT (F6. 0) 
CALL COSTAP (N,X) 


CALL ONE R FT (LOG2N,X) 


DETERMINE THE SPECTRUM ESTIMATES G(P) 


13 


N1=N/2 +1 

J=1 

G (1) =X (1) **2#2.0*0.0001/N 

DO 13 J=2,N1 

K=N+2-J 

G (J) = (X (J) *#24+X (K) **2) *2.0*0 0001/8 
CONTINUE 


SCALE DATA BY 1/0.875 DUE TO COSTAPER 


20 


100 


DO 20 I=1,J3 
G (I) =G (I) /0.875 


CALL FREQ (N,G,J,LB,NR,SAM,N7,VB) 


CONTINUE 

CALL CGPL (XeYoYoNeOe1,1,1,1,HA,HB,HC, VA, VB, VC,ALPH, 6) 
STOP 

END 
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SUBROUTINE COSTAP (N,X) 


REAL X (10) 

N1=N/10 

DO 10 L=1,N1 
A=3.1416/2.0/N1*L 
K=N-L+1 

X (L) =X (L) -X (L) *COS (A) 
X (K) =X (K) ~X (K) *COS (A) 
CONTINUE 

RETURN 

END 


235 


epee GY Sgt Hewat were 


n pee rw 


SES eh, Use 
wee 


or) eee ee « “eu Me | : 


_ oa 
one hee 


a pads We A olen 


236 


SUBROUTINE PREQ(N,G,J,L,NR,SAM,N7, VB) 


wee axe ee ae ere ee eo ee ee ee 


DIMENSION G(J),P(50),AVE(50), ALPHA (20) 
NO=J 
H=0.0001 
NV=0 
DO 20 K=1,NO,L 
NV=NV+1 . 
IF(NV.GT.50) GO TO 49 
A=0.0 
LL=K+L-1 
DO 10 J=K,LL 
10 A=A+G(J) 
AVE(NV)=A/L 
ML=LL-K 
P (NV) = (FLOAT (K)+ (FLOAT (ML) /2.)) /PLOAT (N) /8 
20 CONTINUE 
49 WRITE(6,50) SAM,N7 
50 FORMAT (*1',//,20X,"SPECTRUM ESTIMATES FOR °,A4,//, 
*20X,"RUN NUMBER ',11,//) 
DO 60 JK=1,NR 
60 WRITE(6,55) F(JK),AVE(JK) 
55 FORMAT (15X,F7.0,F10.2) 
READ(5,65) (ALPHA(I) ,I=1,20) 
65 FORMAT (20A4) 
CALL CGPL(F,AVE, F,NRyl51¢ 1,491,020, 100.,8.0,0.0,VB ,6-0,ALPHA,6) 
RETURN 
END 


SUBROUTINE ONE R FT (LOG2N,X) 
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Cocccecccce ee cecece cere cesesecreseerecece eee eseeesesreeevreseeesceeoereee eee 


Cc 
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100 


ONE REAL FOURIER TRANSFORMS 


INTEGER LOG2N 
REAL X (10) 


INTEGER J,JN,K,KN,N,N OVER 2 
REAL ARG,C,PI,S,T, X1,XR,YI,YR 


PI=3.141592654 
N=2** (LOG2N-1) 
CALL R SUB B O 


CALL R SUB B O (LOG2N-1,X(1)) 
CALL R SUB B O (LOG2N-1,X (Nt1)) 
CALL MR 1D FT (LOG2N-1,X(1) ,X(N+1)) 


N OVER 2=N/2+1 


DO 100 J=2,N OVER 2 


K=Nt+2-J 
JN=J+tN 
KN=K+¢N 


XR=(X(J) +X (K)) *.5 
XI=(X (JN) -X(KN))*.5 

YR= (X (JN) +X (KN)) #25 
YI=(X(K) -X (J)) *.5 
ARG=PI*FLOAT (J-1) /FLOAT (N) 


C=COS (ARG) 
S=SIN (ARG) 
T=YR*Ct+YI¥*S 
YI=YI*C-YR*S 
YR=T 

X (J) =XR+YR 

X (K) =XR-YR 

X (KN) =XI+¥I 
X(JN)=YI-XI 
CONTINUE 
XR=X (1) +X (N#1) 
X (1) =XR 

X (N+1) =YR 
RETURN 

END 
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SUBROUTINE RECOVR (LOG2N, RE, IM) 


(Cicraelelelateve elelctelelelevetels lel aiels ele) © 61s 016 (6 0) 6 6) «10 6 a6 6ce) 60 0) 6..010:67 0/610 612: 6)'0'e 6.8) ,0 019 0.0.50. 6,08 


aqnagaagaaanaa 


10 


20 


RECOVERS THE REAL AND THE IMAGINARY PARTS OF THE POURTER 
TRANSFORM WHICH HAS BEEN STORED IN COMPACT HERMETIAN FORM. 


THE COMPACT STORED POIRIER TRANSFORM IS ORIGINALLY IN RE( ) 
SEE TWO R FT FOR DETAILS. 


IN EPFECT THIS SUBROUTINE DOUBLES THE STORAGE AREA REQUIRED 
BUT GIVES EASE IN THE MANIPULATION OF THE COEFFICIENTS WHERE 
FTX(J) = RE(J) + IM(J) 


REAL RE (10), IM (10) 
INTEGER LOG2N, N, N UPON2, MORE 


N = 2**LOG2N 
N UPON 2 = N/2 
MORE = N UPON 2 # 2 


IM(1) = 0.0 

DO 10 J=2,N UPON 2 
I = N-Jt#2 
IM(J) = RE(I) 
CONTINUE 

J = N UPON 2 #1 

IM(J) = 0.0 

DO 20 J=MORE,N 
I = N-J#2 
IM(J) = -RE(J) 
RE(J) = RE(I) 
CONTINUE 

RETURN 


END 
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SUBROUTINE MR 1D FT (LOG2N,X,Y) 
Coccccvcecccc ccc r ee esc ere ese esses eeeecere ee eseceeseereeseeeeoseeeresee eee 
Cc MIXED RADIX ONE DIMENSIONAL FOURIER TRANSFORM 
Cc 
INTEGER LOG2N 
REAL X (10), Y (10) 


INTEGER JJ,J0,31,32,33,N,84 
REALSARG SC lpCes co, l0,l lp he, 2590 RI, RZ yp Roe 5 lp oer5se © 


INTEGER A,B,C,D,E,¥,G,H,1,J,K,L, 4, BS,CS,DS, ES, FS,GS,HS,IS,I3S,KS, 
Toho ieee, Ly FL, PL, Gly hly lL, JL, KL; UL, 9, ot 8)ie) UC 83) 
EQUIVALENCE (BS,S(2)), (CS,S(3)) 4 (DS,S(4)) 4 (ES¢S(5)) » (FSS (6)) 6 
~(GS,S(7)) ¢ (HS¢S (8) ) ¢ (IS¢5 (9) ) ¢ (ISS (10)) > (KS,S(11)) » (LS,S(12)) » 
- (MS,S(13)) ,(AL,U(1)) -(BL,U(2)) ¢ (CL,U(3)) -(DL,U(4)) , (EL, U(5) ) 
-(FL,U(6)) ,(GL,U(7) ) , (HL, U(8)) » (IL, U (9) ) ¢ (JL, U (10)) , (KL, U(11)) 4 

s (LU, G2), (4 70 (13)) 


N=2**LOG2N 
IF (LOG2N.LE.1) GO TO 500 
DO 400 K=2,LOG2N,2 
M=2** (LOG2 N-K) 
N4=4%N 
DO 300 J=1,M 
ARG=6. 28318530 7# FLOAT (J-1) /PLOAT (M9) 
C1=COS (ARG) 
S1=SIN (ARG) 
C2=C1*C1-S1*S1 
S2=C1*S1+C1*S1 
C3=C2*C1-S2*S1 
S3=C2*S14+S52¥*C1 
DO 200 I=M4,N,M4 
JO=I+J3-M4 
J1=J0+s 
J2=J1+8 
J3=J2+M 
RO=X (JO) +X (J2) 
R1=X (J 0) -X (J2) 
IO=¥ (JC) +¥ (32) 
I1=Y (JO) -¥ (32) 
R2=X (31) +X (J3) 
R3=X (J 1) -X (33) 
I2=Y (J 1) +¥ (33) 
I3=Y¥ (J 1) -¥ (J3) 
X (JO) =RO+R2 
Y (J0) =10+12 
IF (ARG.EQ.0.0) GO TO 100 
X (J2) = (R1+13) *C14#(I1-R3) *S1 
Y (J2) = (I1-R3) *C1-(R1+13) #51 
X (J1)= (RO-P2) *C2+ (10-12) *S2 
Y (31) = (10-12) *C2-(RO-R2) *S2 
X (J3) = (R1-13) *C3+ (114+R3) *S3 
Y (J3) = (11+ R3) *C3-(R1-13) *S3 
GO TO 200 

100 CONTINUE 


a 3, a _ : ‘ aa 7 - : l 
oe - ni - PA : : 7 
we ; ; 7 6 ne 
ifs nN 7 ‘. mi 
ae 


or PPrrreuy) Trey ey settee see “s * ” “* e 4 
ty A oe a ee un re 
' 7 : 
‘§ 8 \ y aN aa 
vy on 


2.02 82.8 ast Stata Ouyeresietia 3 ae ‘AA 


.2A,26 e1 2h, @0,eT, ay 20,29 28 Had FG, ee sy Savense: 
{£1} u ers z ft dd, 32, St. Tas 42 un 
« €(2) 2,29) 0010) 25 Bs) (8) 2,80) 4 (EVE 92D) of 
: (UST G2). CUTS, BH), CLOT) R26), (Pp dy | 
eeanaaay (ia uy 36}, ((6)8 49). 10S) 8 ule) s aii) 
(Ct tyo, ii Pr 


Or OT Od (0.0. “08. kia 
rev (Ea-tr verdeteietay ste) iat 
ph ale! baka a- 
Pinice 7 = G)s 

~08%}-SD9 (¢I- tay 
2s (Eger) 45% (E1-7A if : hi ; 
Ee» (E1-th) ED ceaers : 


200 
300 
400 
500 


600 
700 


800 


900 


X (J2)=R1#13 
Y (J2)=11-R3 
X (J1) =RO-R2 
¥(J1)=1I10-12 
X (J3)=R1-13 
¥ (33) =114#R3 


CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 


IF ({LOG2N.EQ.LOG2N/2*2) GO TO 700. 


DO 600 I=1,N,2 
RO=X (I) #X (I+1) 
R1=X (I) -X (I+1) 
IO=¥ (1) +¥ (I+#1) 
I1=¥ (I) -¥ (I+1) 


X (I) =RO 

Y (I) =10 
X (I+ 1) =R1 
Y (I+ 1) =11 
CONTINUE 
CONTINUE 
MS=N/2 
ML=N 


DO 800 K=2,12 


J=14-K 
S(d)c1 


U (J) =S (J+1) 
IF (S(J+1).GT.1) S(J)=S(J+1) /2 


CONTINUE 
AL=BS 


JJ= 


DO 


0 

900 
900 
900 


A=1,AL 

B=A,BL, BS 
C=B,CL,CS 
D=C, DL, DS 
E=D,EL,ES 
F=E,FL, PS 
G=F,GL,GS 
H=G,HL, HS 
I=H,IL,IS 
J=1,0L,0S 
K=J,KL,KS 
L=K, LL, LS 
M=L,ML,MS 


JJ=JJ+1 


IF (JJ.LE.M) GO TO 900 


T=X (JJ) 


X (JJ) =X (M) 


X (M) =T 
T=¥ (JJ) 


¥ (JJ)=¥(M) 


Y(m)=T 
CONTINUE 
RETURN 


240 


BX (Febye= (oye 


one oF OD (8.8 


sh 
' 
oO 


A | on eet : 
» owt die od =p ara) : 
Wa0uTse 


SUBROUTINE R SUB B O (LOG2N,X) 


241 
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INTEGER A,B,C,D,E,¥,G,H,1,J3,K,L,M,N,BS,CS, DS, ES, FS,GS,HS,1S,JS,KS, 


~LS,MS,NS,AL,BL,CL, DL, EL,FL,GL,HL,IL,JL,KL,LL,MNL,NL, 


S (14), U (14) 


POU LVALENCE Bl aspatiz) ), (C5,S (3) )- (Ses (4) ), (85,5 (5) 6 (F5-5 (6) ) 7 
. (65,5 (0er, (85-5 (8) ), (18,519) ), (JS,S (10)) , (KS-5 (11) ) . (L5,5(12)) ¢ 
. (nS, $(13)),(NS,S(14)),(AL,U(1)),(BL,U(2)), (CL,0(3)), (DL,U(4)), 
. (EL, U(5)) »(PL,U(6)), (GL,U(7)) , (HL,U (8) ) , (IL,0(9)) 4 {JL,U(10))e 


- (KL, U(11)) , (LL,U (12) ), (ML,U (13)) - (NL, U (14) ) 


Cc REVERSE SUBSCRIPT BIT ORDER 
Cc 
INTEGER LOG2N 
REAL X (10) 
Cc 
INTEGER JJ 
REAL T 
Cc 
c 
NS=2** (LOG2N-1) 
NL=2*NS 
DO 100 K=2,13 
J=15-K 


U (J) =S (J+1) 


S (J) =1 


IF (S(J+1).GT.1) S(J) =S(J+1) /2 


100 CONTINUE 


AL=BS 


DO 200 
JJ=JJI+1 


IF (JJ.LE.N) GO TO 200 


T=X(JJ) 


A=1,AL 
B=A,BL, BS 
C=B,CL,CS 
D=C,DL, DS 
E=D,EL,ES 
F=E,FL, FS 
G=F,GL,GS 
H=G,HL, HS 
I=H,IL,IS 
J=1I,JL,JS 
K=J,KL,KS 
L=K,LL, LS 
M=L,ML,MS 
N=4,NL,NS 


X (JJ) =X (N) 


X(N) =T 


200 CONTINUE 


RETURN 
END 
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SUBROUTINE ONECAC (LOG2N,X) 
ONE CYCLIC AUTO COVARIENCE 
REAL X (10) 
THIS ROUTINE REPLACES THE POURIER TRANSFORM OF X(STORED IW 
..eCOMPACT HERMITIAN FORM) BY ITS CYCLIC AUTO COVARIANCE 
THE COVARIANCE PUNCTION IS COMPUTED BY TAKING THE INVERSE 
PFT OF PRODUCT OF FFT(X) AND CONJG FFT (X). 


N = 2**LOG2N 
X(1) = X(1) *#2 
NUPON2 = N/2 
STORE ZEROS FOR IMAGINARY PART OF THE PRODUCT. 
DO 100 J = 2,NUPON2 
K = N+2-d 
P = X(J)**2 + X(K) #*2 
X(J) = P 
X(K) = 0.0 
CONTINUE 
K = N/2 + 1 
X(K) = X (K) **2. 


CALL ONECFT (LOG2N, X) 


RETURN 
END 
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SUBROUTINE ONE C FT (LOG2N,X) 
Cee tree ee te nag ale aust eteteiavets’e. 6 ausie aveibie + ees 0 S-aLvie. 07 A) © 9 6 0e0)8 0.0 82 0 Roce aso oie 
G ONE COMPLEX FOURIER TRANSFORM 
c 
INTEGER LOG2N 
REAL X (10) 


INTEGER J,JN,K,KN,N,N OVER 2 
REAL AI,AR,ARG,BI,BR,C,PI,5,T 


N=2** (LOG2N-1) 
PI=3. 141592654 
N OVER 2=N/2+1 
DO 100 J=2,N OVER 2 
K=N+2-J 
JN=J+N 
KN=K+N 
AR=X (J) #X(K) 
AI=X (KN) -X (JN) 
BR=X (J) -X (KR) 
BI=X (JN) +X (KN) 
ARG=PI*FLOAT(J-1) / FLOAT (H) 
C=COS (ARG) 
S=SIN (ARG) 
T=BR¥C+BI*S 
BI=BI*C-BR*S 
BR=T 
X (J) =AR-BI 
X (K) =AR+BI 
X (JN) =BR+ALI 
X (KN) =BR-AI 

100 CONTINUE 

AR=X (1) +X (N#1) 
BR=X (1) -X (#1) 
X (1) =AB 
X (N+1) =BR 
CALL MR 1D FT (LOG2N-1,X(1) ,-X(N*1)) 
CALL R SUB B O (LOG2N-1,X(1)) 
CALL R SUB B O (LOG2N-1,X (N+1)) 
CALL R SUB B O (LOG2N,X) 
RETURN 
END 
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SUBROUTINE PARZEN ( MS, C ) 


Cece eee CRUSOE SH SHSH SEK SHEKHAR HSH SHE CHSC SCE Ke SOSA BES HR SCH SSESO SS OSC 


Cc 
C MODIFIES THE COVARIANCE FUNCTION C(T) BY THE PARZEN LAG WINDOW. 
c THE RESULTS ARE STORED IN C(T), T=1, 2375000 tt? Ve 
Cc 
REAL C (10), D 
INTEGER M, M UPON 2, TWO B 
Cc 


M UPON 2 = M/2 
DO 10 J=1,M UPON 2 
D = 1.0 - 6.0*( FLOAT(J) /FLOAT(M) ) **2 


3 + 6.0*( FLOAT(J)/FLOAT(M) ) **3 
C(J+1) = D * C(J+1) 
10 C Gen T 15N UE 


M UPON 2 = M UPON 2 + 1 
DO 20 J=M UPON 2,M 
D = 2.0*(1.0 - FLOAT(J)/FLOAT(M) ) **3 
C(J+1) = D * C(J+#1) 
20 CO Nott Nios 
RETURYMW 
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SPECIAL COLLECTIONS 


UNIVERSITY OF ALBERTA LIBRARY 


REQUEST FOR DUPLICATION 


I wish a photocopy of the thesis by 


entitled _ 


The copy is for the sole purpose of private scholarly or scientific study 
and research. I will not reproduce, sell or distribute the copy I request, 
and I will not copy any substantial part of it in my own work without per- 
mission of the copyright owner. I understand that the Library performs 
the service of copying at my request, and I assume all copyright responsi- 
bility for the item requested. 


EF. 236 September 1975 


